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Executive  Summary 


In  the  last  few  years  there  has  been  a  renewed  interest  in  the  ultraviolet  (UV) 
spectral  region.  Interest  is  due  to  three  things:  (1)  developments  in  laser 
sources  that  have  considerably  more  power  than  was  available  in  the  past; 
(2)  developments  in  new  detectors  with  excellent  quantum  efficiency,  speed, 
and  the  capability  to  almost  instantaneously  collect  the  full  UV/VIS  spectrum; 
and  (3)  requirements  such  as  remote  sensing  of  biological  agents,  detection  of 
drug  laboratories,  combat  identification  of  friendly  forces,  verification  of 
treaties,  and  detection  of  missile  plumes  and  secure  communications. 

Estimates  of  UV  atmospheric  transmission  are  essential  to  all  requirements  that 
involve  propagation  of  UV  radiation.  UVTRAN  is  an  atmospheric 
transmission  and  laser  radar  (lidar)  model  that  predicts  atmospheric 
transmission  in  the  186  to  700  nm  spectral  region  using  only  visibility,  range, 
^nd  wavelength  as  parameters.  UVTRAN  includes  a  1976  U.S.  standard 
atmosphere  model,  a  modified  Koschmeider  aerosol  model,  and  trace  gases 
(nitrogen  dioxide,  nitrous  oxide,  sulfur  dioxide,  and  water);  the  standard 
atmosphere  ozone  profile  may  be  changed  by  the  user.  The  model  is  designed 
for  any  horizontal,  vertical,  or  slant  path  for  ranges  up  to  11  km  and  for 
visibilities  less  than  50  km.  It  is  primarily  for  continental  air  mass  type  aerosol 
and  is  designed  for  altitudes  up  to  11  km.  A  lidar  model  for  aerosol  elastic 
scattering  or  for  fluorescence  scattering  of  aerosols  or  gases  is  also  part  of  the 
complete  model.  UVTRAN  is  a  Fortran  code  model  for  mainframe  computers 
or  IBM  compatible  PC’s  and  a  Microsoft  Windows  Program  written  in  Visual 
Basic  3.0.  The  Windows  version  has  considerably  more  user  features.  Both 
codes  feature  menu  or  prompt  driven  input  and  plotting  capability. 

UVTRAN  is  an  EOSAEL  92  model  and  is  available  from  the  Army  Research 
Laboratory  to  qualified  requesters.  The  program  is  readily  available  in 
Fortran,  and  the  Visual  Basic  source  code  is  available  in  addition  to  the 
executable  program.  The  Windows  version  comes  complete  with  a  setup 
routine  for  automatic  installation. 
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1.  Introduction 


UVTRAN  is  an  atmospheric  transmission  and  lidar  return  calculation  model  for 
visible  and  ultraviolet  (UV)  wavelengths.  This  EOSAEL  module  combines  a 
transmission  code  suitable  for  use  in  the  visible  and  UV  regions  with  a 
backscatter  code  for  Mie  and  fluorescence  lidar  return  calculations  and  a  sky 
background  radiance  code  into  a  modular,  menu-driven,  user  friendly  program 
for  IBM  PC’s  and  compatible  systems.  The  code  is  available  in  Fortran  source 
code  for  main  frame  computers  and  as  both  source  and  executable  code  for 
IBM  PC’s  and  compatible  systems.  The  PC  version  is  also  available  as  a 
Microsoft  Windows  program  written  in  Visual  Basic.  The  Visual  Basic  version 
has  more  operational  features  and  is  discussed  in  section  4.2. 

The  transmission  model  includes  attenuation  caused  by  molecular  scattering, 
molecular  absorption,  and  particulate  attenuation.  Gaseous  absorption  does  not 
differ  significantly  from  that  in  Lowtran  6;  likewise,  the  calculated  Rayleigh 
scattering  is  only  slightly  less  than  that  of  Lowtran.  The  difference  is  due  to 
matching  of  more  recent  values  of  Rayleigh  attenuation  coefficients.  The 
major  difference  in  the  UVTRAN  and  the  Lowtran  results  are  due  to  the 
different  aerosol  parameterizations  in  the  two  models.  The  Lowtran  models  are 
based  on  aerosol  microphysical  models.  The  UVTRAN  model  is  more  directly 
based  on  empirical  transmission  data,  with  the  wavelength  dependence  of 
aerosol  attenuation  parameterized  in  terms  of  the  visual  range  to  provide  an 
approximate  match  for  visible  and  UV  horizontal  attenuation  data.  The  aerosol 
model  was  compared  with  the  AFGL  standard  aerosol  models  and  with 
experimental  data  on  atmospheric  attenuation  as  a  function  of  the  visual  range. 
Comparisons  are  discussed  in  section  3. 

The  lidar  model  calculates  the  lidar  return  signal  as  a  function  of  range  for  six 
different  lidar  wavelengths.  The  model  is  capable  of  determining  backscatter 
for  Jiominal  conditions  for  Mie  backscatter  caused  by  an  arbitrary  cloud  or 
fluorescent  backscatter  caused  by  particles  or  gaseous  molecules. 
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Default  lidar  system  parameters  are  available,  or  the  user  can  tailor  the  system 
parameters  to  any  particular  system.  Sky  background  is  calculated  for  the 
wavelength  and  set  of  system  parameters  chosen. 


2.  Background 


Electro-optical  sensor  systems  design  requires  a  knowledge  of  the  effects  of  the 
atmosphere  on  the  received  signal  and  on  the  signal  levels  to  be  expected  from 
actively  sensed  atmospheric  targets.  Although  interest  traditionally  focused  on 
infrared  (IR)  and  longer  wavelength  systems,  the  availability  of  improved 
sources  and  detectors  and  the  desire  to  utilize  the  increased  resolution  available 
with  shorter  wavelength  systems  led  to  a  renewed  interest  in  short  wavelength 
optical  systems,  extending  into  the  solar-blind  region  below  290  nm. 

There  has  been  a  need  for  a  simple  code  for  determining  the  parameters  in  the 
visible  and  UV  regions.  To  meet  this  need,  the  UVTRAN  program  was 
developed  to  model  atmospheric  transmission  and  lidar  return  at  visible  and  UV 
wavelengths.  This  model  is  unique  because  it  combines  a  transmission  code 
suitable  for  use  in  the  visible  and  UV  regions  with  a  backscatter  code  for  Mie 
and  fluorescence  lidar  return  calculations  and  a  sky  background  radiance  code 
into  a  modular,  menu-driven,  user  friendly  Fortran  program  for  an  IBM  PC  or 
compatible  system. 

The  transmission  model  differs  from  the  Lowtran  6  [1]  or  similar  models 
because  it  is  designed  for  use  at  visible  and  UV  wavelengths  only,  so  that  the 
band  models  and  emission  calculations  needed  at  IR  wavelengths  are  not 
included.  Molecular  oxygen  absorption  is  included,  a  feature  not  available  on 
PC  versions  of  Lowtran;  absorption  of  trace  gases  nitrogen  dioxide  (NOj), 
nitrous  oxide  (NO),  sulfur  dioxide  (SOj),  and  water  (HjO)  are  also  included. 
In  addition,  the  wavelength  dependence  of  the  aerosol  attenuation  is 
parameterized  in  terms  of  the  visual  range  in  this  model  rather  than  having  a 
wavelength  dependence  independent  of  the  visual  range.  This  different  aerosol 
parameterization  was  chosen,  in  part,  because  of  data  showing  differences 
between  the  Lowtran  results  and  short  wavelength  experimental  data.  [2]  The 
background  radiance  also  is  parameterized  in  terms  of  a  zenith  angle 
dependence  rather  than  being  calculated  on  the  basis  of  single  scattering  theory. 
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This  section  describes  the  computer  program,  discusses  the  scientific  basis  for 
the  model,  and  compares  its  results  with  Lowtran  calculations  and  with 

experimental  data. 

Transmission  Model  Description 

The  transmission  model  determines  the  path  transmission  between  185  and 
700  nm  caused  by  aerosol  attenuation,  molecular  scattering,  and  molecular 
absorption.  The  model  is  designed  to  use  over  relatively  short  ranges  in  the 
troposphere;  consequently,  a  plane  parallel  geometry  oitiy  is  assumed.  The 
model  does  not  include  effects  caused  by  density  gradients  or  atmospheric 
turbulence.  In  all  cases,  the  attenuation  is  assumed  to  be  given  by  the 
Beer-Lambert  attenuation  law: 


where 

a  =  the  attenuation  coefficient 

X  =  the  path  length. 

The  model  does  not  include  a  band  absorption  model  because  of  the  interest  in 
UV  and  visible  absorption  in  which  electronic  absorption  is  important.  No  IR 
water  bands  are  included  in  the  model;  therefore,  the  long  wavelength  limit  for 
the  model  is  approximately  700  nm. 

The  model  is  primarily  designed  for  use  in  the  lower  troposphere. 
Concentrations  of  attenuating  components  are  determined  for  altitudes  of  1 1  km 
or  less.  Because  no  spherical  geometry  is  included,  the  model  is  most 
appropriate  for  relatively  short  ranges. 

The  attenuation  caused  by  molecular  scattering  is  calculated  by  the  molecular 
scattering  formula  based  on  a  fit  to  data  of  Penndorf,  [3]  which  was  used  in 
Lowtran  5  with  a  modification  to  take  account  of  differences  in  the  Lowtran 


formula  and  more  recent  results  in  the  Handbook  of  Geophysics  and  Space 
Environment.  [4] 
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where 

X  =  the  wavelength  of  light 

N  =  the  number  density  of  air  molecules  calculated  from  pressure 

and  temperature  data 

No  =  the  number  density  under  standard  conditions. 

As  written,  the  Penndorf  equation  requires  X  with  units  of  cm  and  calculates 
attenuation  with  units  of  km  *.  UVTRAN  converts  the  entry  units  of  nm  to  the 
wavelength  units  needed  for  the  calculation.  The  correction  term  C, 
C  =  0.987,  is  the  modification  to  the  Penndorf  formula  to  bring  the  calculated 
Rayleigh  scattering  into  agreement  with  the  standard  molecular  scattering  given 
by  the  Handbook  of  Geophysics  and  the  Space  Environment.  [3,4] 

Xlolecular  absorption  in  the  UVTRAN  model  is  calculated  for  the  two  major 
absorbing  species,  oxygen  (O2)  and  ozone  (O3),  and  for  the  minor  absorbing 
species,  NOj,  NO,  SO2,  and  HjO  vapor.  The  choice  of  these  species  for 
inclusion  into  the  model  was  made  on  the  basis  of  expected  concentrations  of 
trace  species  from  the  1976  U.S.  Standard  Atmosphere  [5]  and  molecular 
absorption  coefficients  discussed  in  Ackerman  [6]  and  Calvert  and  Pitts.  [7] 
Under  normal  conditions,  it  is  expected  that  the  molecular  absorption  is  almost 
entirely  determined  by  the  O2  and  the  O3  absorption.  The  other  absorption 
coefficients  are  included  because  high  concentrations  of  these  other  gaseous 
constituents  could  be  present  under  some  circumstances,  leading  to  significant 
absorption. 
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Plots  of  the  absorption  cross  sections  for  O2  and  O3  used  in  UVTRAN  are 
shown  in  figures  1  (O2)  and  2  (O3).  The  model  resolution  is  2.5  nm  for  Ae 
data  and  for  the  O3  Hartley-Huggins  bands.  The  Chappuis  band  data  for  O3 
have  resolutions  of  5  to  10  nm. 


The  resolution  data  is  relatively  low,  but  interpolation  should  not  cause  any 
significant  deviations  from  actual  values,  except  in  the  case  of  the  O2 
Schumann-Runge  bands  between  185  and  200  nm  in  which  the  modeled  values 
are  approximate  averages  over  the  adjacent  bands.  In  this  case,  the  assumption 
was  made  that  measurements  in  this  region  are  relatively  wide  band,  so  the 
averaging  is  appropriate. 


Ozone  Absorption  Data 


Figure  2.  Ozone  absorption  cross  section  in  cm^. 


Aerosol  attenuation  was  parameterized  in  terms  of  visibility  according  to 
equation  (3)  with  the  wavelength  of  light  [8] 

^  [  3.912  * 


550 


(3) 


where 

®^ier(^)  “  desired  particle  extinction 

0^(550)  =  the  Rayleigh  extinction  at  550  nm 
V  =  the  observed  visibility 

q  =  the  Angstrom  exponent  describing  the  wavelength 

dependence  of  the  attenuation. 
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q  is  empirically  related  to  the  visual  range  V  for  V  <  24  km  by 

q  =  .585 

and  for  greater  visual  ranges,  by 

q  =  1.7. 


(4) 


(5) 


Equation  (4)  is  based  on  an  empirical  fit  to  short  wavelength,  relatively  low 
visibility  transmission  data.  [9,10]  TTie  Umiting  case  equation  (5)  is  based  on 
experimental  data  for  high  visibiUty  conditions.  The  relationship  in 
equation  (3)  is  an  extension  of  the  Koschmeider  relation  between  visibility  and 
atmospheric  extinction  at  0.55  /im,  V  =  3.912/a.  [10]  For  V  in  kilometers 
the  units  of  a  are  kilometers  ^ 

Because  the  visibility  depends  on  the  total  atmospheric  extinction  coefficient, 
the  molecular  attenuation  at  550  nm  has  been  subtracted  from  the  total 
atmospheric  attenuation  coefficient  to  determine  the  particle  extinction.  The 
separation  is  made  so  wavelength  scaling  and  particulate  attenuation  variation 
with  height  can  be  made  in  terms  of  aerosol  effects  only. 

Figure  3  shows  a  plot  of  the  calculated  wavelength  dependent  aerosol 
attenuation  for  different  visibiUties.  In  figure  3,  the  dotted  line  corresponds  to 
a  visibility  of  V  =  50  km,  the  solid  line  to  V  =  10  km,  the  dashed  and  dotted 
line  to  V  =  2  km  and  the  dot-dash  line  to  V  =  0.5  km.  The  calculated 
decrease  in  the  wavelength  dependence  of  the  particle  attenuation  with 
decreasing  visibility  is  consistent  with  the  idea  that  the  greatly  reduced 
visibilities  are  associated  with  larger  particles  that  are  the  background  cases  of 

high  visibility. 

Figure  4  shows  a  plot  of  the  total  attenuation  coefficient  as  a  function  of 
wavelength  and  the  attenuation  coefficients  caused  by  aerosol  attenuation, 
molecular  scattering,  and  molecular  absorption  for  assumed  standard 
conditions.  The  total  attenuation  coefficient  is  given  by  the  solid  line,  the 
aerosol  attenuation  by  the  dotted  fine,  the  molecular  scattering  attenuation  by 


the  dashed  and  dotted  line,  and  the  O3  and  Oj  attenuations  by  the  long  and 
short  dashed  lines.  The  calculations  are  made  for  standard  temperature  and 
pressure  conditions  at  sea  level  for  an  assumed  visibility  of  10  km  and  a 
nominal  O3  concentration.  Figure  4  shows  the  relative  importance  of  the 
individual  components. 

In  the  program,  the  pressure  and  temperature  as  a  function  of  altitude  are  given 
by  a  standard  atmosphere  calculation.  Concentrations  of  O3,  SO2,  H2O,  N2O, 
and  NO2  can  be  varied;  O3  altitude  profiles  can  be  entered,  and  altitude 
independent  values  of  the  other  gaseous  constituents  can  be  input.  As 
discussed  above,  the  aerosol  extinction  is  calculated  from  the  visual  range, 
which  is  a  user  input.  The  variation  of  aerosol  extinction  with  altitude  is  done 
in  a  manner  similar  to  that  of  the  Army  vertical  structure  algorithm  in 
Lowtran.  The  percent  change  in  extinction  with  a  given  change  in  altitude  is 
specified  and  used  in  subsequent  calculations. 
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The  transmission  is  calculated  for  one  or  more  ranges  at  a  constant  zenith 
angle.  For  slant  or  vertical  ranges  the  extinction  varies  as  a  function  of  height 
along  the  path;  and  so  the  range  is  divided  into  subintervals,  with  each 
subinterval  covering  an  altitude  range  of  no  more  than  100  m.  The  attenuation 
is  assumed  to  be  constant  in  each  subinterval.  For  multiple  ranges,  the 
transmission  is  done  sequentially  with  the  total  transmission  equal  to  the 
transmission  of  the  final  range  times  the  previous  transmission.  Schematically, 

(6) 
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The  transmission  at  any  level  is  given  by  the  equation 


(7) 


t  = 


e 


-o‘r 


where 


=  the  attenuation  coefficient  at  some  level. 


a  is  the  sum  of  the  attenuation  caused  by  Rayleigh  scattering,  aerosol 
attenuation,  and  molecular  absorption  according  to 


o 


ms 


a  +  o  + 

aer  ozone 


mg 


(8) 


where 


'mg 


=  the  Rayleigh  scattering  of  atmospheric  molecules 
=  the  attenuation  due  to  aerosol  particles 
=  the  absorption  due  to  O3 
=  the  absorption  due  to  the  other  mixed  gases 
=  the  incremental  path  length  through  the  layer. 


2.2  Lidar  Return  Calculations 


The  UVTRAN  module  calculates  the  lidar  return  power  as  a  function  of  range 
for  realistic  lidar  geometries.  Nominal  system  parameters  of  transmitted 
energy,  receiver  area,  receiver  field  of  view  and  spectral  bandwidth,  and 
system  efficiencies  are  provided;  or  alternative  parameters  can  be  input. 
Return  calculations  can  be  made  for  one  of  six  lidar  wavelengths:  532  nm 
(Nd-YAG  doubled),  355  nm  (Nd-YAG  tripled),  266  nm  (Nd-Yag  quadrupled), 
351  nm  (XeF  Excimer),  308  nm  (XeCl  Excimer),  and  248  nm  (KrF  Excimer). 
Calculations  can  be  made  for  elastic  aerosol  and  molecular  backscatter  or  for 
backscatter  cuased  by  fluorescence  from  molecular  or  particulate  material. 
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Elastic  backscatter  is  given  by  the  equation 


P(R)  =  E 


^  4R^ 


where 


=  the  return  power  as  a  function  of  range 
=  the  range 

=  the  Mie  backscatter  coefficient 
=  the  aperture  diameter 
=  the  transmitted  pulse  power 
=  efficiency. 


The  transmission  tx  is  given  by 


-f^AX)rdr 

r,  =  e 


The  fluorescent  backscatter  is  given  by  the  equation 

-  n2  (11) 

mo  = 

with  /Spi,  the  fluorescent  backscatter  coefiicient  for  particulate  fluorescent 
materials  given  by 


Pfl  ” 


where 


=  the  specific  fluorescence  coefficient  of  the  fluorescent 
material  in  the  lidar  spectral  bandwidth 
=  the  mass  concentration  of  the  fluorescing  material, 


18 


and  /Spi,  the  fluorescent  backscatter  coefficient  for  gaseous  fluorescent  materials 
given  by 

.  (13) 

■  4n 

where 

F,,  =  the  fluorescence  cross  section  per  molecule  for  the  spectral 

bandwidth  of  the  lidar 

Cv  =  is  the  concentration  of  the  fluorescing  material  in  ppbv 

p„  =  the  number  concentration  of  air  molecules. 

2.3  Background  Radiance  Calculations 

The  background  radiance  values  can  be  calculated  for  three  cases:  daytime 
clear  sky  illumination,  daytime  cloudy  illumination,  and  nighttime  illumination. 
The  values  are  parameterized  in  terms  of  a  zenith  angle  dependence  relative  to 
assumed  horizon  radiance  levels;  no  solar  azimuth  or  elevation  angle 
dependence  is  assumed.  Sky  and  terrain  backgrounds  are  considered  for 
different  zenith  angles;  sky  for  zenith  angles  less  than  or  equal  to  90®  and 
terrain  for  zenith  angles  greater  than  90®. 

The  background  radiance  routine  radiance  values  for  daylight  conditions  are 
contained  in  a  table  with  radiance  values  at  20  nm  intervals  between  280  and 
700  nm.  Sky  radiance  values  below  280  nm  are  assumed  to  be  0.  The  horizon 
brightness  levels  were  inferred  from  various  data  sources,  including  the 
Infrared  Handbook  and  references  therein.  [11]  The  zenith  angle  dependence 
for  zenith  angles  of  less  than  or  equal  to  90®  for  clear  daylight  conditions  is 
given  by 

W  (019 0-81  sine) 
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and  for  overcast  conditions,  by 


(15) 


W  =  (1-0  2COS0) 


where 

I^r(X)  =  the  nominal  horizon  brightness 

d  =  the  zenith  angle. 


Nighttime  radiance  values  are  assumed  to  be  1.0  x  10'^  times  the  daytime 
values. 


For  zenith  angles  greater  than  90®,  the  brightness  of  the  ground  for  day 
conditions  is  given  by 


4(X)  =  4,r(X)  (1-0  O.7COS0) 


(16) 


The  night  value  is  1.0  x  10*®  times  the  daytime  value. 


3.  Model  Verification 


3.1  Grade  of  Software 

The  atmospheric  transmission  model  should  be  classed  as  developmental 
software;  the  lidar  transmission  calculations  should  be  classed  as  fieldable 

software. 

3.2  Model  Failure 

The  number  of  ranges  called  for  in  the  calculation  should  not  exceed  the 
maximum  number  specified  in  the  interactive  input  file;  otherwise,  dimensions 
would  be  exceeded  and  unpredictable  output  results.  The  maximum  altitude  of 
the  calculation  is  limited  to  11  km  because  flat  earth  geometry  is  assumed  (this 
should  not  be  serious,  the  ranges  for  short  wavelengths  are  gener^y  short). 
A  test  of  mayimiini  altitude  is  included  in  the  program,  which  is  halted  if 
maximum  altitude  is  exceeded.  Other  incorrectly  formatted  input  causes  an 
error  halt  to  the  program. 

3.3  Verification  Tests 

A  comparison  of  our  model  and  the  Lowtran  6  calculations  for  standard 
conditions  is  shown  in  table  1.  The  calculations  were  made  for  wavelengths 
of  300,  550,  and  700  ran,  for  the  atmospheric  conditions  listed  in  the  table. 
The  visibility  was  specified  to  be  10  km,  and  the  wavelength  dependence  of  the 
aerosol  attenuation  is  given  by  equation  (4)  for  UVTRAN  and  by  the  rural 
aerosol  model  for  a  relative  humidity  of  70  percent  for  Lowtran.  [1] 
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Table  1.  Tr.n«n.kdnn  for  UVTRAN  and  Lowtran  for  standard  conditions* 
Wavelength  300  nm  550  nm  700  nm 

T  ^  ^ 

UVT  Low  XJVT  Low  UVT  Low 


r  »> 

.239 

.233 

.892 

.890 

.956 

.957 

•  mn 

.773 

.768 

.998 

.998 

.999 

.999 

‘03 

r 

.013 

.064 

.206 

.206 

.347 

.307 

‘•er 

7-oUi' 

.986' 

.0024 

0.115 

.997* 

.182 

.183 

.331 

.954“ 

.280 

•p  =  1013.25  mb,  T  =  288.15  K,  Oj  =  25.2  ppbv,  Vis  =  23  km,  Range  -  10  km 
=  the  transmittance  caused  by  molecular  scattering 


^To3  =  the  transmittance  caused  by  ozone 

=:  the  transmittance  caused  by  aerosol  particles 
Vo*  =  the  transmittance  caused  by  any  other  gaseous  species 
^caused  by  SO2 


^caused  by  NO2 

^caused  by  water  vapor 

Vjoi  =  the  total  resulting  transmittance 


The  Rayleigh  scattering  attenuation  is  slightly  less  than  the  comparable  Lowtran 
attenuation  as  discussed  above,  leading  to  sUghdy  higher  transmittance  in  the 
model.  Similarly,  there  are  no  significant  differences  in  the  O3  absoiption  in 
the  different  models.  The  calculated  aerosol  attenuation  values  are  significantly 
different  at  300  and  700  nm  in  the  two  models,  with  the  greater  differences  at 
300  nm.  The  agreement  in  the  aerosol  attenuation  at  550  nm  is  due  to  the  fact 
that  the  attenuation  at  this  wavelength  is  determined  by  the  assumed  visual 
range  in  each  of  the  models. 

The  differences  in  the  aerosol  attenuation  in  the  two  models  indicate  that  the 
variation  of  attenuation  with  wavelength  is  less  in  the  AFGL  rural  aerosol 
model  than  in  XJVTRAN.  For  the  cases  shown  above,  q  =  1.7  for  UVTRAN 
compared  with  q  =  0.9  for  Lowtran.  While  the  experimental  data  on  the 
wavelength  dependence  of  the  atmospheric  attenuation  is  Umited,  data  do 
suggest  a  larger  value  of  q  between  300  and  550  nm  than  predicted  by  the 
AFGL  models  except  under  conditions  of  greatly  reduced  visibility. 
Middleton  [10]  reported  that  the  average  value  of  q  was  1.6  under  good  to 
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excellent  seeing  conditions  and  1.3  under  average  seeing  conditions. 
Junge  [12]  used  a  value  of  q  =  1.3  to  describe  the  average  atmospheric 
conditions.  The  experimental  data  of  Wolff,  which  were  used  by  Loehle  [9] 
to  derive  the  empirical  relationship  in  equation  (3),  predict  a  q  less  than  1  but 
are  for  relatively  low  visibility  conditions. 

Baum  and  Dunkelman  [2]  present  a  more  detailed  set  of  wavelength  dependent 
attenuation  data  measuring  attenuation  coefficients  between  250  and  550  nm  in 
Pasadena,  CA  for  a  range  of  visibilities.  They  developed  an  empirical 
relationship  between  visibility  and  wavelength  dependent  atmospheric 
attenuation.  Tables  2  and  3  show  aerosol  attenuation  comparisons  calculated 
by  UVTRAN  and  Lowtran  with  the  attenuation  calculated  from  the  empirical 
relation  of  Baum  and  Dunkelman  [2]  for  visual  ranges  of  300  and  350  nm, 
respectively.  It  is  only  for  visual  ranges  of  less  than  5  km  that  the  q  for 
UVTRAN  is  less  than  the  q  for  the  AFGL  rural  aerosol  model.  The  UVTRAN 
and  Baum  and  Dunkelman  [2]  values  are  closer  at  350  than  300  nm  suggesting 
that  higher  than  standard  O3  concentrations  during  the  Baum  and 
Dunkelman  [2]  measurements  may  be  affecting  the  300  nm  comparison. 


Table  2.  A  comparison  of  observed  and  calculated  attenuation  coefficients  at 
300  nm  for  difierent  visibilities  _ 


Visual 

Range 

(km) 

^tet 

(UVTRAN) 

(kmO 

(Lowtran) 

(km-0 

<^.er 

(B  &  D-emp) 
(km‘) 

100 

.137 

.047 

.213 

40 

.284 

.150 

.340 

20 

.476 

.320 

.553 

10 

.809 

.662 

.978 

5 

1.41 

1.34 

1.83 
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Table  3.  A  comparison  of  observed  and  calculated  attenuation  coefficients  at 
350  nm  for  different  visibilities 

Visual 

Range 

(km) 

Ofex 

(UVTRAN) 

(km‘) 

Ofgf 

(Lx)wtran) 

(km**) 

^tet 

(B  &  D-emp) 
(km**) 

100 

.092 

.042 

.097 

40 

.212 

.134 

.198 

20 

.376 

.288 

.366 

10 

.691 

.594 

.702 

5 

1.21 

1.20 

1.38 

The  comparisons  indicate  that  the  value  of  q  is  related  to  the  visibiUty  and  that 
the  parameterization  for  q  used  in  UVTRAN  is  adequate,  though  certainly  not 
exact  Discussion  of  Woodman  [13]  indicates  that  the  wavelength 
parameterizations  of  equations  (2)  and  (3)  should  be  adequate  throughout  the 
visible  region.  Although  the  relationship  (3)  was  derived  on  the  basis  of  low 
visibility  data  only,  these  comparisons  with  the  Middleton  data  [10]  suggest  the 
model  should  be  adequate  for  high  and  low  visibiUties. 
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4.  Operations  Guide 
4.1  Fortran  Version 


The  MS  DOS  version  of  UVTRAN  is  run  at  the  prompt  by  typing 
UVTRAN  {enter}.  The  program  may  also  be  run  by  typing 
UVTRAN  <  filename. ext  (enter),  filename. ext  is  a  data  file  containing  the 
inputs  to  the  UVTRAN  prompts.  This  latter  method  is  a  property  of  MS  DOS 
that  allows  UVTRAN  to  be  run  from  a  shell  program  in  a  batch  file  and  is 
given  as  an  example  below. 

4.1.1  Inputs 

Input  to  the  program  is  interactive.  Choices  are  presented  for  the 
environmental  parameters  and  for  the  test  conditions  to  be  simulated.  The  user 
must  choose  conditions  or  accept  default  values.  In  general,  letters  may  be 
entered  in  upper  or  lower  case  (i.e.,  when  the  user  is  asked  to  enter  a  choice 
for  a  Mie  or  fluorescence  lidar,  if  the  choice  is  a  Mie  calculation,  m  or  M  may 
be  entered). 

Generally,  entries  are  self  explanatory.  The  aerosol,  elastic  scattering  (Mie) 
lidar  calculation  is  based  on  user  input  of  excess  scattering  material  in  addition 
to  the  scattering  material  estimated  from  the  visibility.  The  excess  material  is 
entered  in  terms  of  backscatter  or  attenuation;  the  nominal  relationship  may  be 
used,  or  a  specific  backscatter  to  attenuation  ratio  may  be  input.  For  the 
fluorescence  lidar,  the  backscatter  is  assumed  to  be  due  to  fluorescence  only. 
Inputs  to  determine  the  fluorescent  backscatter  depend  on  whether  the 
fluorescing  material  is  assumed  to  be  a  gas  or  an  aerosol.  If  the  fluorescent 
material  is  a  gas,  the  fluorescent  cross  section  per  molecule  and  the 
concentration  in  ppbv  must  be  entered.  If  the  fluorescent  material  is 
particulate,  the  specific  fluorescence  cross  section  and  the  concentration  of  the 
fluorescent  material  in  g/m^  must  be  entered. 

Section  5  presents  samples  of  interactive  sessions  for  the  input  of  data  to 
calculate  transmission  or  lidar  return  power  data  shown  in  section  4.2. 
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4.1.2  Output 

Two  sample  output  listings  for  a  transmission  calculation  and  a  LIDAR 
calculation  follow: 

Output  A:  Transmission  Calculation  Output 
UV  transmission  calculations: 


UVTRAN  model  ozone  concentrations  are: 


altitude 

0.  km 

03 

altitude 

= 

1.  km 

03 

altitude 

2.  km 

03 

altitude 

= 

3.  km 

03 

altitude 

= 

4.  km 

03 

altitude 

5.  km 

03 

altitude 

= 

6.  km 

03 

altitude 

= 

7.  km 

03 

altitude 

= 

8.  km 

03 

altitude 

= 

9.  km 

03 

altitude 

= 

10.  km 

03 

concentration 

40.0  ppbv 

concentration 

35.8  ppbv 

concentration 

32.5  ppbv 

concentration 

33.2  ppbv 

concentration 

34.0  ppbv 

concentration 

37.5  ppbv 

concentration 

41.5  ppbv 

concentration 

50.0  ppbv 

concentration 

59.5  ppbv 

concentration 

87.0  ppbv 

concentration 

-  131.0  ppbv 

UVTRAN  model  trace  gas  concentrations  are: 


1  S02  concentration 

2  N20  concentration 

3  H20  concentration 

4  N02  concentration 


1.0  ppbv 
270.0  ppbv 
1000000.0  ppbv 
1.0  ppbv 


visibility  =  10.00  km  at  an  altitude  of  1.00  km 
particle  extinction  changes  by  .0  % 
per  100  meters  altitude  increase 

initial  altitude  =  1.00  km  zenith  angle 


.0  degrees 
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Range 

=  1.00  km 

Final  Altitude 

=  2.00  km 

Wavelength 

Wavelength 

=  500.00  nm 
=  550.00  nm 

Transmission 

Transmission 

=  .641E+00 

=  .676E+00 

Range 

=  2.00  km 

Final  Altitude 

=  3.00  km 

Wavelength 

Wavelength 

=  500.00  nm 
=  550.00  nm 

Transmission 

Transmission 

=  .411E+00 

=  .458E+00 

Output  B:  Lidar  Calculation  Output 
Lidar  return  calculations  for  a  Mie  lidar: 


UVTRAN  modeled  ozone  concentrations  are: 


altitude 

= 

0.  km 

03  concentration 

-  40.0 

ppbv 

altitude 

1.  km 

03  concentration 

-  35.8 

ppbv 

altitude 

= 

2.  km 

03  concentration 

-  32.5 

ppbv 

altitude 

= 

3.  km 

03  concentration 

-  33.2 

ppbv 

altitude 

= 

4.  km 

03  concentration 

-  34.0 

ppbv 

altitude 

= 

5.  km 

03  concentration 

-  37.5 

ppbv 

altitude 

6.  km 

03  concentration 

-  41.5 

ppbv 

altitude 

= 

7.  km 

03  concentration 

-  50.0 

ppbv 

altitude 

= 

8.  km 

03  concentration 

-  59.5 

ppbv 

altitude 

= 

9.  km 

03  concentration 

-  87.0 

ppbv 

altitude 

= 

10.  km 

03  concentration 

-  131.0 

ppbv 

UVTRAN  modeled  trace 

1  S02  concentration 

2  N20  concentration 

3  H20  concentration 

4  N02  concentration 


concentrations  are: 

1.0  ppbv 
270.0  ppbv 
1000000.0  ppbv 
1.0  ppbv 


gas 


visibility  =  10.00  km  at  an  altitude  of  1.00  km 


particle  extinction  changes  by  1.0 

per  100  meters  altitude  increase 


Lidar  system  parameters: 

Pulse  Energy  =  1.00  joules; 

Receiver  fov  =  2.00  mrad; 

transmit  efficiency  =  .55; 


Receiver  Diameter 
Spectral  Width 
receiver  efficiency 


.600  m; 
.50  nm; 
.54 


initial  altitude 


=  1.00  km  zenith  angle 


=  30.0  d^rees 


added  values: 

range  #  value  added  for  beta 


1  .200E-05 


Lidar  wavelength  = 

532.0  nm 

bkg  signal  = 

.264E-07  watts 

range 

[km] 

final  alt 
[km] 

beta(tot) 

[(msr)**l] 

beta  ratio 

2xtrans 

sogtot 

[watts] 

sigtKay; 

[watts] 

1.00 

2.00 

1.87 

2.73 

.180E-04 

.159E-04 

.134E-I-00 

.129E-I-02 

.426E-I-00 

.168E+05 

.966E-04 

.843E-05 

.722E-05 

.654E-06 

4.2  Visual  Basic  for  Windows  Version 

The  purpose  of  the  Windows  version  of  UVTRAN  is  to  give  the  end  user  a 
user  friendly  graphical  interface  and  provide  several  operation^  features  not 
available  through  Fortran  such  as  online  help  and  built-in  plotting  capability. 
This  has  been  done  by  porting  the  Fortran  source  code  over  to  Microsoft 
Visual  Basic.  UVTRAN  for  Windows  is  a  complete  Windows  based 
application  designed  and  written  within  the  Microsoft  Windows  environment. 
XJVTRAN  for  Windows  requires  Microsoft  Windows  version  3.0  or  later  to 
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For  specific  hardware  and  operating  system  requirements  see  the  Microsoft 
Windows  User’s  Guide.  A  Visual  Basic  runtime  module,  VBRUN100.DLL 
(or  VBRUN300.DLL  depending  on  the  latest  version  under  which  UVTRAN 
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for  windows  was  compiled)  must  be  present  for  it  to  execute  outside  of  Visual 
Basic.  The  module  is  supplied  with  UVTRAN  for  Windows. 

4.2.1  Installation 

The  Windows  version  of  UVTRAN  comes  with  a  setup  program  on  the 
diskette.  The  user  inserts  the  diskette  into  the  disk  drive;  from  Program 
Manager  in  Windows,  puts  the  cursor  on  File;  and  from  the  pull  down  menu, 
clicks  on  Run.  In  the  pop-up  box,  the  user  types  a:setup  {enter}.  The 
program  has  defaults  for  installation  and  tells  how  much  hard  disk  space  is 
required.  The  default  directory  may  be  changed.  The  installation  program 
creates  a  group  with  the  appropriate  icons  and  the  appropriate  directory  on  the 
hard  disk  and  decompresses  the  files  to  the  directory. 

4.2.2  Initiating  a  Session 

From  Program  Manager  in  Windows,  the  user  opens  the  group  containing  the 
UVTRAN  for  Windows  icon  and  executes  the  program  by  selecting  the 
UVTRAN  for  Windows  icon.  As  with  all  Microsoft  Windows  applications, 
UVTRAN  for  Windows  can  be  minimized  or  quit  anytime  the  control  box  is 
available  to  the  user.  For  a  complete  discussion  of  suspending  and  stopping 
the  Windows  application,  please  refer  to  the  Microsoft  Windows  user  guide. 

UVTRAN  for  Windows  was  designed  to  minimize  data  input  errors  and  make 
the  program  easier  to  use.  The  major  convention  used  is  the  placement  of 
input  buttons  and  data  display.  To  enter  new  data  the  user  simply  clicks  on  the 
proper  button  and  an  input  window  will  open  for  the  user  to  input  the  new 
value.  If  data  values  are  present  in  the  current  window  and  buttons  are  visible, 
the  user  clicks  on  the  data  item  to  be  updated.  Each  window  containing  data 
that  the  user  can  edit  or  update  contains  an  Accept  Values  button.  This  button 
saves  the  values  display  in  the  window  and  returns  the  user  to  the  previous 
window.  Because  of  the  design  of  UVTRAN  for  Windows,  there  is  no  need  for 
ordered  procedures  to  run  the  model.  All  system  parameters  are  set  to  default 
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values  at  start  up,  allowing  the  user  greater  flexibility  to  process  data  using  a 
what  if  methodology.  A  discussion  of  each  window  and  the  data  displayed  in 

them  follows: 

UVTRAN  for  Windows  Opening  window.  This  is  the  main  menu/window  for 
the  UVTRAN  for  Windows  model.  This  window  displays  five  buttons  for  the 
user  to  choose  from:  (1)  Run  it  (or  Process),  (2)  Setup  Aerosol  Profiles, 
(3)  Edit/Set  Default  Ozone  Levels,  (4)  Edit/Set  Default  Trace  Gas  Levels,  and 
(5)  Calculation  Type.  The  Run  it  (or  Process)  button  runs  the  model  using  the 
system  parameters  that  the  user  sets  in  the  other  four  windows  accessed 
through  the  buttons  available  on  the  Opening  window. 

Setup  Aerosol  Profiles.  Selecting  this  button  displays  the  Aerosol  Profiles 
window.  The  Aerosol  Profiles  window  displays  the  parameters  pertaining  to 
the  aerosol  profile  used  in  the  model.  The  user  can  change  and  edit  the  aerosol 
profile  used  by  selecting  different  actions  from  this  window.  When  saUsfied 
with  the  selected  value,  the  user  cUcks  on  Accept  Values  to  close  the  window 
and  return  to  the  previous  window. 

Edit/Set  Default  Ozone  Levels.  Selecting  this  button  selects  the  Ozone  ^vels 
window.  The  user  can  edit  any  of  the  O3  levels.  When  satisfied  with  the 
selected  value,  the  user  cHcks  on  Accept  Values  to  close  the  window  and  return 
to  the  previous  window.  The  user  may  also  click  on  the  Reset  Default  Values 
button  to  return  the  O3  levels  to  their  default  values. 

Edit/Set  Default  Trace  Gas  Levels.  Selecting  this  button  displays  the  Trace 
Gas  Concentrations  window.  The  user  can  edit  any  of  the  trace  gas  values. 
When  satisfied  with  the  selected  value  the  user  cUcks  on  Accept  Values  to  close 
the  window  and  return  to  the  previous  window.  The  user  inay  also  click  on 
the  Reset  Default  Values  button  to  return  the  O3  levels  to  their  default  values. 

Calculation  Type.  Selecting  this  button  displays  the  Calculation  T^  window. 
In  this  window  the  user  selects  Transmission  or  Lidar  calculation  from  the 
model.  The  user  sets  the  number  of  ranges  for  the  model  to  compute  as  well 
as  the  shortest  distance  and  length  of  interval  of  measurements.  The  user  may 
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edit  system  parameters  for  the  select  calculation  type  or  accept  default  values. 
If  Edit  Parameters  is  selected,  the  proper  window  appears  allowing  the  user  to 
set  parameter  values  for  Transmission  or  Lidar  calculations. 

Transmission  Setup.  If  Transmission  was  selected  in  the  Calculation  Type 
window,  the  user  will  see  the  Transmission  Setup  displayed  when  Edit 
Parameters  is  selected  from  the  Calculation  Type  window.  The  user  may 
set  the  number  of  wavelengths  from  this  window  and  enter  specific  values 
for  the  desired  wavelengths. 

lidar  Setup.  If  Lidar  was  selected  in  the  Calculation  Type  window,  the 
Lidar  Setup  is  displayed  when  the  user  selects  Edit  Parameters  from  the 
Calculation  Type  window.  The  user  may  set  all  lidar  parameters  from  this 
window,  including  special  values  for  Mie  or  Fluorescence  type  lidar 
calculations. 

Run  It!  (or  Process).  Selecting  this  button  executes  the  model  using  the 
parameters  set  in  the  setup  windows. 

Results.  When  all  calculations  are  complete  the  Results  window  is  displayed 
showing  the  user  the  resulting  calculation  from  the  input  parameters  passed  to 
the  model  and  giving  the  user  the  choice  to  print  the  results  or  place  them  in 

a  file. 
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5.  Sample  Data  Input 


The  following  exEmples  show  an  interactive  session  for  transmission  and  Lidar 
calculations.  Computer  output  queries  and  information  are  shown  in  normal 
typeface;  user  responses  are  shown  in  bold  type,  and  comments  on  the  response 
are  shown  in  italics.  In  some  cases,  changes  are  entered  that  are  the  same  as 
the  default  values  to  illustrate  procedures. 


The  following  input  generates  the  Fortran  version  sample  run  for  the 
transmission  calculations  shown  in  section  4.  User  responses  are  in  bold  and 
program  writes  to  the  screen  are  in  italics. 

Do  you  want  hard  copy  output?  (y  or  n)  y 

Make  sure  printer  is  on. 

present  values  for  trace  constituents  are  based 
on  1976  standard  atmosphere  values.  Ozone  data 
are  given  for  each  km  for  altitudes  between  0 
and  20  km.  Other  trace  gas  concentrations  are 
representative  lower  tropospheric  values.  All 
values  are  given  as  number  mixing  ratios. 

Enter  return  for  default  parameters;  o  to 
change  ozone  parameters,  t  to  change  any  of  the 
other  trace  gas  concentrations;  or  a  to  change 
both  ozone  and  other  trace  gases.  {return} 

The  drfault  visual  range  is  10  km 
Hit  s  for  d^ault  value 

else  enter  the  desired  visual  range  (in  km)  s 
enter  altitude  of  receiver  (in  km)  1.0 
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enter  altitude  of  visibility  measurement 
enter  s  if  same  as  receiver  altitude  s 

Visual  range  selected  is  10  km 
at  an  altitude  of  1 .0  km 

enter  number  corresponding  to  desired  aerosol  profile 

1  if  aerosol  attenuation  is  constant  with  height 

2  if  aerosol  attenuation  increases  with  height 

3  if  aerosol  attenuation  decreases  with  height 

(2  or  3  assume  a  change  by  a  factor  of  x  %/y  m— 

a  change  of  -0.97%/100  m  corresponds  to  an  approximately  constant  mass 
mixing  ratio)  1 

Is  path  horizontal  (h),  vertical  (v),  or  slant  (s)  v 
enter  u  or  d  for  upward  or  downward  looking  systems  u 
Do  you  want  transmission  or  LIDAR  calculation  (Enter  Tor  L)  T 
How  many  ranges  for  calculation  (  <  =  100  )  2 


enter  shortest  range  (km)  1.0 

enter  range  interval  (km)  1.0 

Should  transmission  be  high  (h)  or  low  (1)  resolution 
(h  only  for  laser  transmission)  1 

Enter  number  of  wavelengths  (  <  20)  for  transmission  calculations  2 
enter  wavelength  values  in  nm  500.  550. 

{The  Program  then  calculates  the  transmission  for  the  selected  wavelengths  and 
ranges.  The  remainder  of  output  is  an  echo  of  the  ou^ut  to  the  printer  above.} 


The  following  Fortran  version  sample  run  will  generate  the  sample  run  for  the 
lidar  calculations  shown  in  section  4: 

Do  you  want  hard  copy  output?  (y  or  n)  y 

Make  sure  primer  is  on. 

present  values  for  trace  constituents  are  based 
on  1976  standard  atmosphere  values.  Ozone  data 
are  given  for  each  km  for  altitudes  between  0 
and  20  km.  Other  trace  gas  concentrations  are 
representative  lower  tropospheric  values.  All 
values  are  given  as  number  mixing  ratios. 

Enter  return  for  default  parameters;  o  to 
change  ozone  parameters,  t  to  change  any  of  the 
other  trace  gas  concentrations;  or  a  to  change 
both  ozone  and  other  trace  gases,  [return] 

The  d^ault  visual  range  is  10  km 
Hit  s  for  default  value 

else  enter  the  desired  visual  range  (in  km)  s 

enter  altitude  of  receiver  (in  km)  1.0 

enter  altitude  of  visibility  measurement 
enter  s  if  same  as  receiver  altitude  s 

Visual  range  selected  is  10  km 
at  an  altitude  of  1.0  km 

enter  number  corresponding  to  desired  aerosol  profile 

1  if  aerosol  attenuation  is  constant  with  height 

2  if  aerosol  attenuation  increases  with  height 

3  if  aerosol  attenuation  decreases  with  height 

(2  or  3  assume  a  change  by  a  factor  of  x  %/y  m~ 
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a  change  of -0.97% /1 00  m  corresponds  to  an  approximately  constant  mass 
mixing  ratio)  2 

enter  percent  change  in  attenuation  and  vertical  distance  (in  meters) 
over  which  change  occurs  (positive  numbers)  1.0  100. 

Is  path  horizontal  (h),  vertical  (v),  or  slant  (s)  s 

enter  zenith  angle  for  slant  path  30.0 

Do  you  want  transmission  or  UDAR  calculation  (Enter  Tor  L)  L 

How  many  ranges  for  calculation  (<  -  100)  2 

enter  shortest  range  (km)  1.0 

enter  range  interval  (km)  1.0 

Enter  background  condition  for  calculation 
daylight  clear  =  1 

daylight  overcast  =  2 

night  =31 

Clear  Daylight  Background  Selected 

Udar  return  calculation  can 

be  made  for  6  different  laser  wavelengths 

1  ~  532  nm  (Nd-YAG  doubled) 

2  —  355  nm  (Nd-YAG  tripled) 

3  —  266  nm  (Nd-YAG  quadrupled) 

4  -  351  nm  (XeF  Excimer) 

5  -  308  nm  pCeCl  Excimer) 

6  -  248  nm  (KrF  Excimer) 

Enter  the  number  of  the  desired  wavelength  1 

Number  chosen  is  1  wavelength  chosen  is  532.0  nm 
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Following  parameters  are  default  parameters  for  Lidar  system: 


1.  Laser  Pulse  Energy  --  1.00  joules 

2.  Receiver  Mirror  dicaneter  --  0.6  m 

3.  Receiver  Field  of  View  —  2. 00  mrad 

4.  Spectral  Bandwidth  of  System  —  0.5  nm 

5.  Transmitter  Efficiency  —  .55 

6.  Receiver  Efficiency  —  -54 

Enter  number  of  parameter  to  change; 

or  enter  s  if  current  values  are  used.  1 


Enter  new  value  of  parameter  2.00 

Do  you  want  to  change  another  parameter  (y  or  n)  y 

Current  Udar  parameters  are 

1.  Laser  Pulse  Energy  —  2.00 Joules 

2.  Receiver  Mirror  diameter  ~  0.6  m 

3.  Receiver  Field  of  View  ~  2.00  mrad 

4.  Spectral  Bemdwidth  of  System  —  0.5  nm 

5.  Transmitter  Efficiency  —  ‘55 

6.  Receiver  Efficiency  —  -54 

Enter  number  of  parameter  to  change; 

or  enter  s  if  current  values  are  used.  1 

Enter  new  value  of  parameter  1.00 

Do  you  want  to  change  another  parameter  (y  or  n)  it 
Do  you  want  Mie  or  fluorescence  calculations?  (M  or  F)  M 
Mie  lidar  was  selected 

Program  assumes  a  beta  profile(m**-l  sr**-l)  determined 
from  aerosol  and  molecular  profiles.  An  additional  component 
can  be  added  at  any  point  in  the  profile  by  inserting  beta 
or  extinction  values  for  each  range. 


Enter  number  of  ranges  at  which  excess  aerosol  is  to  be  added  1 
Do  you  want  to  enter  as  backscatter  or  attenuation  (b  or  a)?  b 
Do  you  want  to  enter  a  specific  beta/attenuation  ratio  (y  or  n)  n 
Ranges  and  numbers  listed  below. 

Enter  range  number  and  value  (ordered  pair)  for  each  point 
at  which  excess  aerosol  is  added.  Attn=pcm-1]  beta— [(msr)-l] 

#  range  it  range  tt  range  tt  range 

1  1.000  2  2.000 

1  2.0e-6 

The  printed  output  is  also  output  to  the  screen.  The  program  then  calculates 
the  lidar  return  for  the  selected  ranges.  The  remainder  of  output  is  an  echo  of 
the  output  to  the  printer  above. 


6.  Conclusions 


UVTRAN  is  an  effective  model  for  quick  estimates  of  UV  transmission  and 
lidar  returns.  The  algorithms  used  to  determine  the  transmission  and  the  lidar 
returns  have  been  verified,  and  UVTRAN  compares  favorably  with  the  Baum 
and  Dunldeman  data  and  model.  [2]  It  also  compares  favorably  with 
Lowtran  6  when  compatible  scenarios  are  selected.  The  model  is  not  intended 
to  replace  more  complete  transmission  models  such  as  Lowtran  and  full-blown 
radiative  transfer  models;  however  it  can  estimate  transmission  for  situations 
described  in  this  report. 


UVTRAN  is  a  very  easy  model  to  use  and  has  a  very  friendly  user  interface- 
especially  the  Windows  version.  The  Windows  version  also  comes  with  an 
automatic  setup  installation  program.  This  model  is  available  through  the 
EOSAEL  library  from  ARL  or  from  the  authors.  The  Fortran  source  code  is 
included,  and  the  Visual  Basic  3.0  code  is  available  upon  request. 
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Appendix 

UVTRAN  Fortran  Listing 


c 

c 

c 

c 

c 

c 

c 

c 

c 


Program  UVTRAN 

Version  2.1  June  18,  1994 

Program  calculates  uv  and  visible  wavelength  transmission, 
and  LIDAR  return.  Program  requires  as  input  the  wavelength(s), 
visibility,  range,  zenith  angle,  altitude  above  sea  level, 
concentration  of  gaseous  absorbers.  Program  will  output 
either  transmission  or  expected  LIDAR  return  together  with 
an  estimate  of  the  background. 


common/ilas/ilas 

common/za/zref,alfapt,alfsgn,att55r 

common/o3con/ o3con(2 1)  ,o3pmt(2 1) 
common/amolcn/  so2con  ,an2ocn ,  h2ocon ,  ano2cn ,  trppb(4) 
common/zl/zl(21) 
common/syspar/syspar(8) 

^  j|e  3K  3(e  *  JK  **’*'*  ♦*’•'***  * 

c  identification  of  character  variables 

c  trflg— transmission  or  lidar  calculation 

c  pthtype-’horizontal,  vertical,  or  slant  path 
c  udflag-up  or  down  looking  system 

c  resflg — high  or  low  resolution  transmission 

c  prton-'print  or  not  flag 

c  pflg-change  any  lidar  system  parameters 

c  tyflg--Mie  or  fluorescence  lidar 

c  pgflg--particle  or  gas  for  fluorescence 

c  dflag— change  ozone,  other  trace,  or  all 

c  iflg2--flag  to  review  trace  gas  concentrations 

c  iflg3--flag  to  change  trace  gas  concentrations 

c  baflg-flag  to  indicate  wheter  beta  or  attn  is  input 
c  bsflg--flag  to  indicate  whether  bseta/sigma  ratio  is  imput 

Jle**  ♦  **•«**♦**♦*****  **^**  **  ** 

character  trflg,  pthtype,  udflag,resflg,prton 
character  pflg,  tyflg,  pgflg,  baflg,  bsflg 
characterise  msga(9) 
character  dflag,iflg2,iflg3 
character*56  msgb(7) 
parameter  (nrmax = 100) 


dimension  alamb(20),trans(20,nrmax) 
dimension  hrlamb(6) 
dimension  trcon(4) 

dimension  range  (nrmax),  zf(nrmax),betrat(nrmax) 
dimension  betadd(nrmax),  attadd(nrmax),betfl (nrmax) 
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dimension  betray(nnnax),  betpar(nrmax),  betsum(nrmax) 

dimension  sigsum(nrmax),  sigray(nrmax),  sigfl(nrmax) 

dimension  tr2x(nrmax) 

dimension  flcon(nrmax) 

dimension  imum(nrmax),  addval(nrmax) 

e<}uivalence(trcon(  1 ) ,  so2con) ,  (trcon(2)  ,an2ocn) ,  (trcon(3)  ,h2ocon) 

&,(trcon(4),ano2cn),(trppb(l),so2ppb),(trppb(2),an2opb) 

&,  (trppb(3)  ,h2oppb) ,  (trppb(4)  ,ano2pb) 


data  cdtrad/.017453292/ 

data  hrlamb/532. ,355. ,266. ,35 1 . ,308. ,248./ 

data  vis/ 10.0/ 

data  msga/ 

r  High  resolution  transmission  calculations  can 
be  made  for  6  different  laser  wavelengths 


T 

y 

4’ 

5’ 

6’ 

T 

8’ 
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1  -  532  nm  (Nd-Yag  doubled) 

2  —  355  nm  (Nd-Yag  tripled) 

3  -  266  nm  (Nd-Yag  quadrupled) 

4  -  351  nm  (XeF  Excimer) 

5  -  308  nm  (XeCl  Excimer) 

6  —  248  nm  (KrF  Excimer) 

9’  Enter  the  number  of  the  desired  wavelength 

data  msgb/ 

r  enter  number  corresponding  to  desired  aerosol  profile  ’, 

2’  1  if  aerosol  attenuation  is  constant  with  height  ’, 

3’  2  if  aerosol  attenuation  increases  with  height 

4’  3  if  aerosol  attenuation  decreases  with  height  ’, 

5’  (2  or  3  assume  a  change  by  a  factor  of  x  %/y  m- 

6’  a  change  of  -0.97%/ 100m  corresponds  to  an  approximately’, 

T  constant  mass  mixing  ratio.  ’/ 


c 

c  Initialize  and  clear  screen 
c 

call  clear 

write(*,’(////a//a///a///)’) 

&’  PROGRAM  UVTRAN’, 

&’  Version  2.1’, 

&’  Program  calculates  transmission  or  Lidar  returns.’ 

c 

c 

c  Input  if  hard  copy  output  is  wanted, 
c 

write(*,’(///a\)’)’  Do  you  want  hard  copy  output?  (y  or  n) 
read(*,’(al)’)prton 
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ipmt=0 

if(prton.eq.’y’ .or.prton.eq.  ’  Y’)then 

ipmt=l 

w^ite(^’(///a)’)’  ’ 

pause’  Make  sure  printer  is  on.  ’ 

w 

c  Geophysical  and  environmental  set  up 
c 

c  initialization  routine  for  gaseous  constituents 

c  subroutine  checks  to  see  if  default  values  of  ozone  and  other 
c  trace  gases  are  ok,  if  not,  the  user  is  given  the  opportunity 
c  to  change  the  parameters  in  tcondata. 
c 

c  set  up  ppb  values  for  printing 
c 

do  630  i=  1,21 
o3pmt(i) =o3con(i)*  1  .Oe+9 
630  continue 
do  634  i=  1,4 
trppb(i)  =  trcon(i)*1.0e+9 
634  continue 

c  determine  whether  mods  are  needed  in  o3  or  other  concentrations 

write(*,’(/a/)’)’  ’  •  *  k  c»h’ 

write(*,’(a)’)’  present  values  for  trace  constituents  are  based  ^ 

write(*!’(a)’)’  on  1976  standard  atmosphere  values.  Ozone  data 
write(*/(a)’)’  are  given  for  each  km  for  altitudes  between  0 
write(*’,’(a)’)’  and  20  km.  Other  trace  gas  concentrations  ^e’ 
write(*,’(a)’)’  representative  lower  tropospheric  values.  ^  All 
write(*,’(a)’)’  values  are  given  as  number  mixing  ratios  , 
write(*!’(/a)’)’  Enter  return  for  default  parameters;  o  to  ’ 
write(*,’(a)’)’  change  ozone  parameters,  t  to  change  any  of  me 
write(*,’(a)’)’  other  trace  gas  concentrations;  or  a  to  change 
write(’^!’(a)’)’  both  ozone  and  other  trace  gases’ 

650  read(*,’(al)’)dflag 
if(dflag.eq.’  ’)  go  to  500 

if(dflag.eq.’o’.or.dflag.eq.’0’)go  to  660 

if(dflag.eq.’t’.or.dflag.eq.’T’)go  to  700 
if(dflag.eq.’a’.or.dflag.eq.’A’)go  to  660 
write(*,’(a\)’)’  Please  hit  return  or  enter  o,  t,  or  a 
goto  650 
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c 

c 

c 


Make  ozone  modifications 


660  write(*,’(/a/)’)’  current  model  ozone  concentrations  are: 

write(*  ,20 1 0)(z  1  (i)  ,o3pmt(i)  ,i =1,21) 

665  write(*,’(/a)’)’  enter  altitude  (0,1,2,. ..20  km)  and  new  ozone  con 

Icentration  (ppbv)’ 

670  read(’'',*)zin,o3inpt 

if  (zin.ge.0.0.and.zin.le.20.)  go  to  680 
write(*,’(/a'a\)’)’  e>^ler  alt.  in  km  increments  between  0  and  20’, 
r  reenter  alt  and  new  ozone  concentration’ 
goto  670 

680  j=  int(zin+0.5) 

o3con(j  + 1)  =o3inpt’^  1  .Oe-09 
write(*,2015)j  ,o3inpt 

write(*,’(a\)’)’  do  you  want  to  review  data  set?’ 
read(*,’(a)’)iflg2 

if(iflg2.eq.  ’y’  .or.iflg2.eq.  ’  Y’)then 
do  674  i=  1,21 
o3pmt(i) = o3con(i)  1 .  Oe + 9 
674  continue 

write(*,’(/a/)’)’  current  model  ozone  concentrations  are:’ 
write(*,  1010)(zl(i),o3pmt(i),i= 1,21) 
endif 

write(*,’(a\)’)’  do  you  want  to  change  another  ozone  value?’ 

685  read(V(a)’)iflg2 

if(iflg2.eq.’y’.or.iflg2.eq.’Y’)  goto  665 
if(iflg2.eq.’n’.or.iflg2.eq.’N’)  goto  690 
write(*,’(a\)’)’  need  to  enter  either  y  or  n’ 
goto  685 

690  if(dflag.eq.’o’.or.dflag.eq.’0’)  go  to  500 


up  for  modification  of  trace  gas  concentrations 
c 

700  write(*,2020)(trppb(i),i=l,4) 

730  write(*,’(a)’)’  If  change  in  assumed  concentration  is  desired’ 

735  write(*,’(aH’  enter  the  number  of  trace  component  and  the’ 
write(*,’(a)’)’  new  concentration.  Enter  negative  id  and  dummy’ 
write(’^,’(aH’  concentration  if  no  changes  are  desired’ 
read(* ,  ♦)idno,trppb(idno) 
if(idno.lt.0)go  to  5()0 
trcon(idno)  =trppb(idno)*  1  .Oe-9 
write(*,2025)idno,trppb(idno) 
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write(*,’(a\)’)’  do  you  want  to  review  data  set?’ 

read(*,’(a)’)iflg2 

if(iflg2 .  eq.  ’  y  ’  .or.  iflg2 .  eq.  ’  Y’)then 
write(*,1020)(trppb(i),i=  1,4) 

endif  .  ,  «  , 

write(*,’(a\)’)’  do  you  want  to  change  another  value/ 

750  read(*,’(a)’)iflg3 

if(iflg3.eq.’n’.or.iflg3.eq.’N’)  go  to  500 
if(iflg3.eq.’y’.or.iflg3.eq.’Y’)  go  to  735 
write(*,’(a\)’)’  please  enter  y  or  n’ 
goto  750 

2010  formate  altitude  =  ’,f3.0,’  km  03  concentration  =’,f7.ie  ppbV) 
2015  formate  new  ozone  data:’/’  alt  =  ’,i3,’  km  03  value  ,fo.l, 

1’  ppbv’) 

2020  fomiat(’  current  modeled  trace  gas  concentrations  are.  / 
r  1--S02  concentration  ppbv’// 

2’  2--N20  concentration  ’,f7.1,’  ppbv’// 

3’  3--H20  concentration  ppbv’// 

4’  4--N02  concentration  ’,f7.1,’  ppbv’) 

2025  format(’  ID  for  trace  gas  is  ’,i2,’  cone  =  ’,f7.2,’  ppbv’) 

^:tc*itcm****************************************** 

c  This  section  sets  up  the  profile  of  aerosol  extinction. 

c  Three  choices  are  offered:  . .  .  •  u 

c  1.  Aerosol  attenuation  is  assumed  to  be  constant  with  height 

c  in  the  region  of  interest.  •  u  •  u 

c  2.  Aerosol  attenuation  is  assumed  to  increase  with  height 

c  by  a  factor  of  x%  per  y  m.  .  •  u 

c  3.  Aerosol  attenuation  is  assumed  to  decrease  with  height 
c  by  a  factor  of  X%  per  y  m.  A  change  of  0.97%/100m 

c  will  produce  an  aerosol  profile  with  a  roughly  constant 

c  mixing  ratio. 

c  An  exponential  profile  is  assumed.  The  profile  exponent  is 
c  given  by  alfapt  times  delta  z  (delz). 


500  continue  .  ,a  ,  , 

write(*,’(//a/a/a\)’)’  The  default  visual  range  is  10  km 

&’  Hit  s  for  default  value;’, 

&’  else  enter  desired  visual  range  (in  km). 
read(*,*,  err =502)  vis 
go  to  503 

502  vis=  10.0 

503  continue 
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write(*,’(///a\)’)’  enter  altitude  of  receiver  (in  km)  ’ 
read(’^,’(f6.0)’)  zO 

write(*,’(//a/a\)’)’  enter  altitude  of  visibility  measurement 
1  enter  s  if  same  as  receiver  altitude 
readC"^,  *,err =5 12)zref 
go  to  513 

512  zref=zO 

513  write(*,3020)vis,zref 
call  stdatm(zref,t,p,den) 

att55r  =  (3.912/vis)-den*4.505*1.0e-22 

write(*, ’(//7(a/))’)(msgb(i),i  =  1 ,7) 

520  read(*,’(i3)’)iaprof 
if(iaproheq.  1)  goto  550 
if(iaprof.eq.2)  goto  570 
if(iaprof.eq.3)  goto  570 
write(*,’(a\)’)’  enter  either  1,2,  or  3’ 
go  to  520 

550  alfapt=0.0 
alfsgn=0.0 
go  to  103 

570  write(*,’(/a/a\)’)’  enter  %  change  in  attenuation  and  vertical  dis 
Itance  (in  meters)’,’  over  which  change  occurs  (positive  numbers)’ 
read(*,*)attpct,zdelt 
alfapt=(alog(1.0+0.01*attpct))/zdelt 
if(iaprof.eq.2)  alfsgn=  1.0 
if  (iaprof.eq.3)  alfsgn=  -1.0 
go  to  103 

103  pctlOO  =  (exp(alfsgn*alfapt*100.)  -  l.)*100. 

3020  format(/’  Visual  range  selected  is  ’,f6.2,’  km’,/, 
r  at  an  ^titude  of  ’,f6.2,’  km’) 

^:^**i^*****#m*********************************** 


write(*,’(/7a\)’)’  Is  path  horizontal  (h),  vertical  (v),  or  slant 
&  (s)?  ’ 

105  read(*,’(al)’)  pthtype 

if  (pthtype.eq.’h’.or.pthtype.eq.’H’)  goto  120 
if  (pthtype.eq.’v’.or.pthtype.eq.’V’)  goto  140 
if  (pthtype.eq.’s’.or.pthtype.eq.’S’)  goto  160 
write(*,’(//a\)’)’  Pathtype  not  identified  properly.  Enter 
&  s,  V,  or  h’ 
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goto  105 
c 

c  set  up  for  horizontal  path 
c 

120  thetaz  =  90.0 

thetzr  =  thetaz*cdtrad 
go  to  170 
c 

c  set  up  for  vertical  path 
c 

140  write(*,’(//a\)’)’  Enter  u  or  d  for  upward  or  downward  looking 

&systems  ’ 
read(*,’(bn,al)’)udflag 
thetaz  =  0.0 

if  (udflag.eq.’d’. or. udflag.eq.’D’)  thetaz  =  180. 

thetzr  =  thetaz^cdtrad 
go  to  170 
c 

c  set  up  for  slant  path 

c  , 

160  write(*,’(a\)’)’  enter  zenith  angle  for  slant  path 

read  (*,*)thetaz 
write(*,1007)thetaz 
1007  format(’  zenith  angle  is  ’,f7.2) 
thetzr = thetaz*cdtrad 


170  continue 


^^c**c**m*********************************** 


write(*,’(//a)’)’  ’ 

write(*,’(a\)’)’  Do  you 


want  transmission  or  LIDAR  calculation? 


&(Enter  T  or  L)’ 

200  read(*,’(bn,al)’)trflg 

if  (trflg.eq.’t’.or.trflg.eq.’T’.or.trflg.eq.’l’.or.trflg.eq.  L ) 


&go  to  210 

write(*,’(a\)’)’  a  "t"  or  an  "1"  should  be  entered 


goto  200 


210  write(*,’(//a\)’)’  How  many  ranges  for  calculation  (<  or  -  100)’ 
read(*,*)nrange 
if(nrange.eq.  1)  then 
write(*,’(//A\)’)’  enter  range  (km) 
read(*,*)range(l) 
else 

write(*,’(//A\)’)’  enter  shortest  range  (km) 
read(*,*)range(l) 

write(*,’(//A\)’)’  enter  range  interval  (km) 
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read(*,*)delmg 

endif 

zf(l) =zO+range(l)*cos(thetzr) 
do  230  j  =2,nrange 
rangeO)  =range(  1) +0-  l)*delmg 
zf(j)=zO+range(j)*cos(thetzr) 

230  continue 

if  (trflg.eq.’t’.or.trflg.eq.T’)  go  to  220 
if  (trflg.eq.’l’.or.trflg.eq.’L’)  go  to  400 

220  write(*,’(//a,/a\)’)’  Should  transmission  be  high(h)  or  low(l)  re 
&solution’,’  (h  only  for  laser  transmission)  ’ 
read(*, ’(bn,al)’)resflg 
ires=-l 

if(resflg.eq.  ’1’  .or.resflg.eq.  ’L’)ires=0 
if(resflg.eq. ’h’ .or.resflg.eq. ’H’)ires = 1 
if(ires.ne.0.and.ires.ne.l)go  to  220 

if(ires.eq.0)  then 

write(*,’(////a\)’)’  Enter  number  of  wavelengths  (<  20)  for  trans 
Amission  calculations  ’ 
read(*,*)  nlamb 

write(*,’(///a\)’)’  enter  wavelength  values  in  nm  ’ 
read(*,  *)(alamb(i),i = 1  ,nlamb) 
c 

c  calculate  transmission  values  for  each  of  the  lambda  values 
c  for  each  of  the  ranges 
c 

do  280  i  =  1,  nlamb 

if(alamb(i).lt.  185.0)  alamb(i)=  185.0 
if(alamb(i).gt.700.)  alamb(i)=  700.0 
do  240  j  =  1,  nrange 
if  (j-eq.l)  then 

call  trmod(alamb(i),zO,vis,range(l),thetzr,t,p,den, 

,  trans(i,l),ires) 
else 

call  trmod(alamb(i),zf(j-l),vis,delmg,thetzr,t,p,den, 

,  trans(i,j),ires) 

trans(i,j)  =  trans(i,j)  ♦trans(i,j-l) 
endif 

240  continue 
280  continue 


else  . 

write(*,’(//2(a/),/6(a/)a\)’)(msga(i), 1  =  1,9) 

read(*,*)li 

write(*,  1020)li,hrlamb(li) 

c  combine  the  resolution  parameter  with  the  laser  id  parameter  so  that 
c  ires=o  for  low  resolution  case  and  ires  =  1,2,...  for  the  individual 

c  laser  cases 
ires=li 
nlamb=l 

do  290  j  =  1,  nrange 
if  (j.eq.l)then 

call  tr mod(hrlamb  (li) ,  zO ,  vis ,  range(  1 ) ,  thetzr ,  t ,  p , 

,  den,trans(l,l),ires) 
else 

calltrmod(hrlambai),zf(j-l),vis,delmg,thetzr,t,p,den, 

,  trans(l,j),ires) 

trans(l,j)  =  trans(l,j)  *  trans(l,j-l) 
endif 

290  continue 
endif 

c  write  transmission  data  to  printer 
c 

If(ipmt.eq.  1)  then 
310  continue 

open(4,file=  ’pm’) 

c  open(4,file  =  ’  ’, status  =’new’) 

write(4,’(//a//)’)’  uv  transmission  calculations  ’ 
write(4,’(/a/)’)’  current  model  ozone  concentrations  are:’ 
write(4’,  1008)(zl(i),o3pmt(i),i= 1, 1 1) 
write(4, 1009)(trppb(i),i = 1 ,4) 

write(4,1010)  vis,zref,pctl00 

if  (ires.ne.O)  write  (4,’(a/)’)’  high  resolution  data’ 

write(4,1013)  zO,thetaz 

Do  315  jfl,  nrange 

write  (4,1015)  range(j),  zf(j) 

if(ires.ne.O)  then 

write(4, 1014)hrlambai),trans(l  ,j) 
else 

write(4,1014)(  alamb(i),trans(i,j),  i=l,nlamb) 
endif 

315  continue 
endif 
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c 


write(*,1010)  vis,zref,pctlOO 

if  (ires.ne.O)  write  (*,’(a/)’)’  high  resolution  data’ 

write(*,1013)  zO,thetaz 

Do  318  j  =  l,  nrange 

write  (*,1015)  rangeO),  zf(j) 

if(ires.ne.O)  then 

write(*,  1014)hrlanib(li),trans(l,j) 
else 

write(*,1014)(  alamb(i),trans(i,j),  i=l,nlanib) 
endif 

318  continue 

1008  format(’  altitude  =  ’,f3.0,’  km  03  concentration  - 
1  ppbv’) 

1009  format(/’  current  modeled  trace  gas  concentrations  are.  // 
r  1--S02  concentration  ’,f7.1,’  ppbv’/ 

2’  2--N20  concentration  ’,f7.1,’  ppbv’/ 

3’  3--H20  concentration  ’,f9.1,’  ppbv’/ 

4’  4--N02  concentration  ’,f7.1,’  ppbv’) 

1010  format(/’  visibility  =  ’,f6.2,’  km  at  an  altitude  of  \ 
&f5.2,’  km’/’  particle  extinction  changes  by  ’,f6.1,’  %’ 

&/’  per  100  meters  altitude  increase’) 

1013  format(//lh  ,’  initial  altitude  =  ’, 

&,f5.2,’  km  zenith  angle  =  ’,f5.1,’  degrees’//) 

1014  format(lh  ,’  Wavelength  =  ’,f7.2,’  nm  Transmission  =  ’ 

&,e8.3)  ^  , 

1015  format  (Ih  ,/,’  Range  =  ’,  f5.2,’  km  Final  AlUtude  =  , 

&f5.2,’  km’/) 

1020  format(’  Number  chosen  is  ’,i2,’  wavelength  chosen  is  ’, 
&f5.1,’  nm’) 


c  do  LIDAR  c^culations 
c 

^l^*SII2IIlI®****»****^***»».-**>^*>^»*’^-^********************** 

c  determine  background  condition:  1=  clear  daylight;  2=  overcast 

c  3=  night 

c 

write(*,’(/a/a/a/a\)’)’  Enter  background  condition  for  calculaUo 
&n’,’  daylight  clear  =  1’,’  daylight  overcast  =  2  , 
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&'  night  =  3 

410  read(*,*)  ibkflg 

if(ibkflg.eq.l.or.ibkflg.eq.2.or.ibkflg.eq.3)  goto  412 
write(*,’(//a\)’)’  Please  enter  1,  2,  or  3 
goto  410 

412  if(ibkflg.eq.l)  then 

write(*,’(//a/)’)’  Clear  Daylight  Background  Selected 
endif 

if(ibkflg.eq.2)  then 

write(*,’(//a/)’)’  Overcast  Daylight  Background  Selected 
endif 

if(ibkflg.eq.3)  then 

write(*,’(//a/)’)’  Night  Background  Selected’ 

c 

c  determine  wavelength  of  laser  source 


write(*,’(//a/)’)’  Lidar  return  calculations  can’ 
write(’^’,’(a/,/6(a/)a\)’)(msga(i),  i=2,9) 


read(’*‘,’^)li 

c  determine  whether  default  lidar  parameters  are  ok  as  is  or  whether 
c  changes  should  be  made 


c 


write(*,’(//a//)’)’  Following  are  default  parameters  for  Lidar  sy 


&stem:  ’ 

write(*,1041)(syspar(i),i=l,6)  ^ 

1041  formate  1.  Laser  Pulse  energy  -- ’,f4.2,’ joule /, 

2.  Receiver  Mirror  Diameter  --  ’,f4.3,’  m  ’/ 

3.  Receiver  field  of  view  --  ’,f4.2,’  mrad’/, 

4.  Spectral  Bandwidth  of  System  -  ’,f4.2,’  nm  ’/, 

5.  Transmitter  Efficiency  — 

6.  Receiver  Efficiency  —  ’,f3.2,/, 

Enter  number  of  parameter  to  change;  ’/, 
or  enter  s  if  current  values  are  used.’/) 


2’ 

3’ 

4’ 

5’ 

6’ 

T 

8’ 


404  read(*/,err=406)isp 

write(*,’(/7a\)’)’  Enter  new  value  of  parameter 
read(*,’^)syspar(isp) 

write(*,’(//a\)’)’  Do  you  want  to  change  another  parameter  (y  or 
&n)  ’  ’ 

read(*,’(al)’)pflg 
if(pflg.eq.’y’.or.pflg.eq.’Y’)  then 
write(’^,’(//a/)’)’  Current  Lidar  parameters  are  ’ 
write(*,  1041)(syspar(i),i= 1,6) 
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go  to  404 
else 

L/ 

c  set  up  Lidar  eqn  parameters 
c 

406  continue 

enpls=syspar(l) 

colar=  (syspar(2)’'‘*2)’'‘(3. 141593/4.0) 

recfov  =  3.141593*(0.001*syspar(3)/2)**2 

specbw=  syspar(4) 

efftot  =  syspar(5)*syspar(6) 

eqna  =  enpls*colar*efftot 

iflg=-l 

write(*,’(//a\)’)’  Do  you  want  Mie  or  fluorescence  calculations? 
&(MorF)  ’ 

408  read(*,’(al)’)tyflg 

if(tyflg.eq.’m’.or.tyflg.eq.’M’)  iflg=l 
if(tyflg.eq.’f  .or.tyflg.eq.’F’)  iflg=0 
if(iflg.eq.-l)  then 

write(*,’(//a\)’)’  Please  enter  either  M  orF  ’ 
go  to  408 
else 
endif 


j^jg  lidar,  assume  that  the  backscatter  is  in  three  parts, 
^nc»»ii«*iie|p^;^]g/^ii1ar,  nominal  particle,  and  added  as  an  additional  aerosol  cloud. 

fQr  the  molecular  is  calculated  from  atmospheric  density,  beta 
nominal  particle  is  calculated  from  visibility,  additional  aerosol 
entered  in  terms  of  beta  and  extinction  calculated  for  each. 
c******Nominally,  beta  =  0.03  sigma**.? 

j.#****#**#******#*************************************************’"’"**** 

if(iflg.eq.  1)  then 

write(*,’(//a/a/a/a/a//a\)’)’  Mie  lidar  was  selected.’, 

&’  Program  assumes  a  beta  profile  (m**-l  sr**-l)  determined  ’, 

&’  from  aerosol  and  molecular  profiles.  An  additional  component’, 

&’  can  be  added  at  any  point  in  the  profile  by  inserting  beta’, 

&’  or  extinction  values  for  each  range.’, 

&’  Enter  number  of  ranges  at  which  excess  aerosol  is  to  be  added  ’ 
read  (*,*)  nadd 
if  (nadd.ne.O)  then 

write  (*,’(/a\)’)  ’  Do  you  want  to  enter  as  backscatter  or  attenua 
&tion  (b  or  a)?  ’ 

420  read(*,’(al)’  )baflg 
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if(baflg.eq.’b’.or.baflg.eq.’B’)  i3flg  =  0 
if(baflg.eq.’a’.or.baflg.eq.’A’)  i3flg  =  1 
if(i3flg.eq.-l)  then 

write(*,’(//a\)’)’  Please  enter  an  a  or  a  b. 
go  to  420 

write(*,’(//a\)’)’  Do  you  want  to  enter  a  specific  beta/atten 

, ratio?  (  y  or  n)’ 

421  read(*,’(al)’)bsflg 
i4flg  =  -1 

if(bsflg.eq.’n’.or.bsflg.eq.’N’)  i4flg  -  0 
if(bsflg.eq.’y’.or.bsflg.eq.’Y’)  i4flg  =  1 
if(i4flg.eq.-l)  then 

write(*,’(//a\)’)’  Please  enter  either  y  or  n. 
go  to  421 
endif 

if  (i4flg.eq.  1)  then 

write(*,’(//a\)’)’  Enter  multiplier  for  sigma(km**-l)  to  get  beta 

,  ((msr)-l)  ’ 
read  (*,*)bsrat 
endif 

write  (*, ’(/a/a/a/a/)’)’  Ranges  and  numbers  Usted  below.’, 

&’  Enter  range  number  and  value  (ordered  pair)  for  each  point , 

&’  at  which  excess  aerosol  is  added.  Attn=[km-1]  beta=[msr-l]  , 
&’  #  range  #  range  #  range  #  range  ’ 

write  (*,’(i6,f9.3,i6,f9.3,i6,f9.3,i6,f9.3) ) 

&  (  j,rangeO),  j  =  l,  nrange) 
do  424  i=l,nadd 
if  (iSflg.eq.O)  then 
read  (*,*)  m,  betadd(m) 
imum(i)=m 
addval(i) =betadd(m) 
if  (i4flg.eq.l)  then 
attadd(m)  =  (l/bsrat)*betadd(m) 
else 

attadd(m)  =  148*(betadd(m))**1.43 
attadd(m)  =  attadd(m)* 1000.0 
endif 
else 

read(*,*)  m,  attadd(m) 
imum(i)  =  m 
addval(i)  =  attadd(m) 
if  (i4flg.eq.l)  then 
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betadd(m)  =  bsraf^attadd(m) 
else 

betadd(m)  =  0.03*(attadd(m))**0.7 
betadd(m)  =  betadd(m)*0.001 
endif 
endif 

424  continue 
endif 

do  428  j  =  1,  nrange 

r^ll  betamod(hrlamb(li),zO, vis, ranged), thetzr,t,p, den, 

&betrayO),  betpard),li) 
betsumd)  =  betrayd)  +  betpard)  +  betaddd) 
betratd)  =  betsumd)/betrayd) 

428  continue 

call  trmod(hrlamb(li),zO,vis,range(l),thetzr,t,p,den,tral ,li) 

tr2x(l)  =  trainral 
do  432  j  =2,  nrange 

call  trmod(hrlanib(li),zfd-l),vis,delmg,thetzr,t,p,den, 

&tral,li) 

if  (attaddd-l).ne.O.O)  then 
deladd  =  attaddd-l)*delmg 
tradd  =  exp(-deladd) 
tral  =  tral  *  tradd 
endif 

tra2  =  traPtral 
tr2xd)  =  tra2*tr2xd-l) 

432  continue 


c  complete  Lidar  return  calculations 
c 

do  436  j  =  1, nrange 

eqnb  =(eqna*tr2xd)/((ranged)*1000.)*’"2))*(1.5e+8) 

sigsumd)  =  (eqnb*betsumd)) 
sigrayd)  =  (eqnb*betrayd)) 

436  continue 

call  bkgmd(thetzr,hrlamb(li),ibkflg,bkgval) 
sigbkg  =  bkgval*colar*recfov*specbw 

else 


fluorescence  lidar  assume  that  the  backscatter  is  due  to 
fluorescence  only,  but  that  the  attenuation  is  due  to 


molecular  and  nominal  aerosol. 

l^^^^^^^^^J!***:.*******************************^^^ 

write(*,’(//a/a/a/a/a//aya\)’)’  Fluorescence  lidar  was  selected’,^ 
&’  Program  assumes  an  extinction  profile  (m**-l)  determined 
&’  from  aerosol  and  molecular  profiles.  Beta  is  due  to  ,  ^ 

&’  fluorescence  only,  which  must  be  input.  Fluorescence  can  , 


&’  be  added  at  any  point  in  the  profile.’, 

&’  Do  you  want  to  model  gas  or  particle  fluorescence?  , 

&’  Enter  g  or  p  : 


440  read(*,’(al)’)pgflg 

i2flg  =  -1  .  n 

if(pgflg.eq.’g’.or.pgflg.eq.’G’)  i2flg  -  0 

if(pgflg.eq.’p’.or.pgflg.eq.’P’)  i2flg  =  1 
if(i2flg.eq.-l)  then 

write(*,’(//a\)’)’  Please  enter  a  g  or  a  p. 
go  to  440 
else 
endif 

write(*,’(//a\)’)’  What  is  the  wavelength  shift  of  the  return?  ( 
&in  nm)  ’ 
read(*,*)flshft 
alam2  =  hrlamb(li) + flshft 


if(i2flg.eq.O)  then 
write(*,’(//a/)’) 

&’  What  is  the  cross  section  per  nm  (m**2/nm) 
else 


for  fluorescence? 


write(*,’(//a/)’) 

&’  What  is  the  specific  cross  section  per  nm  (m*’''2/g/nm)  for  iluor 
Aescence?’ 
endif 


read(*,*)  fixe 
flxc2 =flxc*syspar(4) 

write  (*,’(/a\y)  ’  Enter  number  of  ranges  at  which  fluorescence 
&  material  concentration  is  input 
read  (*,*)  nadd 

do  442  j  =  l,nrange 
flconC)  =  0.0 
442  continue 

write  (♦, ’(/a/a/a/a/a//)’) 

&’  Ranges  and  numbers  are  listed  below.’, 

&’  Enter  range  number  and  concentration  value  (as  ordered  pair)  , 
&’  for  each  point.  Enter  as  ppbv  for  gas  concentration  or  ’, 

&’  g/m**3  for  particle  concentration.’. 
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&’  #  range  #  range  #  range  U  range 

write  (*,’(16, f9.3,i6,f9.3,i6,f9.3,i6,f9.3)’) 

&  (  j,rangeO),  j  =  l,  nrange) 

do  444  i=l,nadd 
read  (*,*)  m,  flcon(m) 
imum(i)  =  m 
addval(i)  =  flcon(m) 

444  continue 


caU  trmod(hrlamb(li),zO,vis,range(l),thetzr,t,p,den,tral ,li) 
call  trmod(alam2,z0,vis,range(l),thetzr,t,p,den,tra2,0) 
tr2x(l)  =  tral*tra2 
do  450  j  =  2,  nrange 

call  trmod(hrlamb(li),zf(j-l),vis,delmg,thetzr,t,p,den, 

&tral,li) 

call  trmod(alam2,zf(j-l),vis,delmg,thetzr,t,p,den,tra2,0) 
tral2  =  tral*tra2 
tr2x(j)  =  tral2*tr2x(j-l) 

450  continue 
c 

fluorescettcc  Lidar  return  calculations 
c 

do  454  j  =  1,  nrange 
if(i2flg.eq.O)  then 
call  stdatm(zf(j),t,p,den) 
flcon2  =  flcon(j)*den*1.0e-3 
betfl(j) =flxc2*flcon2/(4*3. 141593) 
else 

betfl(j)  =  flxc2*flcon(j)/(4*3. 14159) 
endif 

eqnb  =eqna*tr2x0)/((range0)**2)*(1.5e+2)) 
sigfl(j)  =  equb^betflCj) 

454  continue 

call  bkgmd(thetzr,alam2,ibkflg,bkgval) 
sigbkg  =  bkgval*colar*recfov*specbw 

endif 

c  write(*,1042)  flsxc2,flmcon,flshfl 

cl042  format(//’  fluorescence  cross  section  (in  wavelength  interval) 
c  &  ’,el0.3,’m**2/g7’  concentration  =  ’,el0.3,’  g/m**3’ 
c  &/’  wavelength  shift  =  ’,f5.1,’  nmV) 


c  writeC^,1043)  flxc2,flppb,flshft 

cl043  format(//’  fluorescence  cross  section  (in  wavelength  interval)  - 
c  &  ’,el0.3,’m**27’  concentration  =  ’,el0.3,’  ppbv’ 
c  &/’  wavelength  shift  =  nmV) 


V' 

If(ipmt.eq.l)  then 
continue 

c  open  (4,  file=’pm’) 

open  (4, file  =  ’  status  =’new’) 
write(4,’(//a\)’)’  Lidar  return  calculations’ 
if(iflg.eq.l)  then 
write(4,’(a/)’)’  for  a  Mie  Lidar’ 
else 

write(4,’(a/)’)’  for  a  fluorescence  Lidar’ 
endif 


write(4,’(/a/)’)’  current  modeled  ozone  concentrations  are;’ 
write(4!  1008)(zl(i),o3pmt(i),i= 1, 1 1) 
write(4,  l(X)9)(trppb(i),i=  1 ,4) 

write(4, 1010)vis,zref,pctl(X) 
write(4, 1046)(syspar(i),i=  1,6) 
write(4,l()47)  zO,thetaz 
if  (iflg.eq.l)  then 
if  (nadd.eq.O)then 

write  (4, ’(a)’)’  no  added  attenuation  values’ 
else 

write(4,’(a,/,a\)’)’  added  values:’,’  range  #  value  added  for’ 
if  (Bflg.eq.O)  then 
write  (4, ’(a)’)’  beta’ 
else 

write  (4, ’(a)’)’  attenuation’ 

endif  .  .  , 

write(4,’(i5,5x,ell.3)’)(imum(i),addval(i),  i  -  l,nadd) 

endif 

write(4,1048)  hrlamb(li),  sigbkg 
write  (4, ’(a, a/a, a/)’)’  range  final  alt  beta(tot)  ’, 

&’  beta  ratio  2xtrans  sigtot  sig(Ray)’, 

&’  [km]  [km]  [(msr)**-l] 

&’  [watts]  [watts]  ’ 

Do456j  =  l,nrange 

write(4, 1044)  range0),zf0),betsum0),betrat0),tr2x0), 
&sigsumO),  sigrayCi) 

456  continue 
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else 


if  (nadd.eq.O)then 

write  (4, ’(a)’)’  no  added  attenuation  values’ 
else 

write(4,’(a,/,a\)’)’  added  values:’,’  range  value  added  for 
.fluorescence’ 

write(4,’(i4,6x,ell.3)’)(imum(i),addval(i),  i  =  l.nadd) 

endif 

write(4,1049)  hrlamb(U),  alam2,  sigbkg 
write(4,’(a)’)’  range  final  alt  beta  (fl)  2xtrans  ret 

&um  signal  ’ 

write(4,’(a/)’)’  [km]  [km]  [(msr)-l] 

&[watts]  ’ 
do  458  j  =  l.nrange 

write(4, 1052)  range(j),zfO),betflC),tr2xO),sigfl(j) 

458  continue 
endif 

else 

endif 

1044  format(lx,f5.2,4x,f5.2,2x,5(el2.3)) 

1045  formatj//’  Lidar  transmitter  wavelength  =  ’,f5.1,’  nm  ’,/ 

&’  fluorescence  return  wavelength  =  ’,f5.1,’  nm’/ 

&’  fluorescence  beta  in  bandwidth  =  ’,e8.3,’  m**-l  sr**-l  / 

&’  total  transmission  =’,e8.3,///, 

&’  return  fluorescence  signal  =  ’,e8.3,’  watts’,/ 

&’  background  signal  =  ’,e8.3,’  watts’) 

1046  formatf//’  Lidar  system  parameters:’/’  Pulse  Energy  = ’, 

&  f4.2,’  joules;  Receiver  Diameter  =  ’,f5.3,’  m;’/ 

Receiver  fov  =  ’,f4,2,’  mrad;  Spectral  Width  =  ’,f4.2, 

&’  nm;’/’  transmit  efficiency  =  ’,f4.2,’;  receiver  efficiency  = 
&’,f4.2/) 

1047  format(/lh,’  initial  altitude  =  ’,f5.2, 

&’  km  zenith  angle  =  ’,f5.1,’  degrees’/) 

1048  format  (/’  Lidar  wavelength  =  ’,f5.1,’  nm  bkg  signal  =  ’ 

,ell.3’  watts’/)  , 

1049  format  (//’  source  wavelength  ’,f5.1,’  nm  return  wavelength  , 

,f5.1’  nm’/’  background  signal  =  ’,ell.3/) 

1052  format(lx,f5.2,4x,f5.2,2x,4(el3.3)) 


c 

c 

c 


if(iflg.eq.l)  then 
endif 
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write(* ,  ’  (//a\)  ’)  ’  LID AR  return  calculations’ 
if(iflg.eq.l)  then 

write(*,’(a/)’)’  for  a  Mie  LIDAR’ 
else 

write(*,’(a/)’)’  for  a  fluorescence  LIDAR’ 
endif 

write(*,’(/a/)’)’  current  modeled  ozone  concentrations  are:’ 
write(*!  1008)(zl(i),o3pmt(i),i=  1 . 1 1) 
write(*,  1009)(trppb(i),i=  1 ,4) 

writeC^,  1010)vis,zref,pctl(X) 
write(*,1047)  zO,thetaz 

if(iflg.eq.l)  then 

write  (*,’(a,a/a,a/)’)’  range  final  alt  beta(tot)  , 

&’  beta  ratio  2xtrans  sigtot  sig(Ray)’, 

&’  [km]  [km]  [(msr)**-l] 

&’  [watts]  [watts]  ’ 

Do  466  j  =  l,nrange 

write(*,  1044)  range0),zf(j),betsum0),betrat0),tr2x0), 
&sigsum(j),  sigrayO) 

466  continue 
else 

write(*,1049)  hrlamb(li),  alam2,  sigbkg 
write(*,’(a)’)’  range  final  alt  beta  (fl)  2xtrans  ret 

&um  signal  ’ 

writeC^/Ca/)’)’  [km]  [km]  [(msr)-l] 

&[watts]  ’ 

do  468  j  =  l,nrange 

write(*,  1052)  range0),zf0),betfl(j),tr2x0),sigfl0) 

468  continue 
endif 


end 

c 

subroutine  trmod(alamb,zO, vis, range, thetzr,t,p, den, trans, ires) 

c  subroutine  trmod  is  the  transmission  module  that  calculates  transmission 
c  over  a  given  range.  Input  parameters  are 
c  alamb-wavelength 
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c  zO  —  receiver  altitude 

c  viS“Visual  range 

c  range— distance  to  target 

c  thetzr-zenith  angle  in  radians 

c  t,p, den-atmospheric  temperature,  pressure,  and  density 

c  trans— returned  transmission 

c  ires-flag  for  resolution  and  for  laser  type;  ires=0  low; 

c  ires=l,  2,  3,  ...  high-  for  laser  line  designated  by  numbers 

c 

nmgnt  =  int((range/0.1)*abs(cos(thetzr)))+l 
mgint  =  (range/nmgnt) 

sigmg  =  0.0 
attmg  =  0.0 
do  275  j  =  1,  nmgnt 

zi  =  zO  +((mgint/2)+((j-l)*rngint))*cos(thetzr) 
call  stdatm(zi,t,p,den) 
call  o3ccal(zi) 
if(ires.eq.O)then 

call  absmol(alamb,den,amolab) 
else 

call  hrabsm(ires,den,amolab) 
endif 

call  attmol(alamb,den,amolsc) 
call  attaer(zi,alamb,den,vis,paitex) 
c  write(*,  ’(3el0.4)’)amolab,amolsc,partex 
sigmg  =  amolab+amolsc+partex 
attmg  =  attmg  +  sigmg*mgint 
275  continue 

trans  =  exp(-attmg) 

return 

end 

-##i|e***#**#*****»**#>|c*********»**********’'‘*’''********’'“''*******’''**’^***’'‘******* 

c 

subroutine  betamod(alamb,zO, vis, range, thetzr,t,p, den, 

&betray,  betpar,li) 
c 

c  subroutine  calculates  beta  values  for  normal  extinction  and 
c  molecular  profiles.  Beta  values  are  in  (m  sr)**-l 
c 

zf  =  zO  +  range*cos(thetzr) 
call  stdatm(zf,t,p,den) 
call  attmol(alamb,den,amolsc) 
betray =(amolsc/8. 130)*!  .Oe-3 
call  attaer(zi,alamb,den,vis,partex) 
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betpar  =  (0.03*partex’*‘*0.7)»1.0e-3 

return 

end 

subroutine  clear 
c 


character  clr 
clr=char(27) 
write(*,100)  clr 
100  format  (lx,al,’[2J’) 
return 


subroutine  bkgmd(thetzr , alamb ,  ibkflg ,  bkg val) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Purpose  is  to  calculate  background  in  watts/ (sr  nm  m**2) 
thetzr  =  zenith  angle  in  radians 
ibkflg  =  flag  which  indicates  conditions  of  illumination 

1  =  daylight  clear  sky 

2  =  daylight  overcast 

3  =  night 


Array  Hbkg  contains  nominal  daylight  horizon  brightness 
values  from  280  to  720  nm.  Otheie  values  are  calcualted 
from  these  values 


Kt************************************** 


dimension  hbkg(23) 

data  hbkg/1.0e-8,4.0e-3,2.3e-2,3.9e-2,5.6e-2,6.5e-2,7.5e-2 
&8.6e-2,9.5e-2,1.01e-l,1.04e-l,1.04e-l,1.01e-l,9.9e-2,9.8e-2, 
&9.2e-2,8.8e-2,8.1e-2,7.3e-2,6.5e-2,4.8e-2,3.6e-2,2.6e-2/ 


if(alamb.lt.280)  then 

bkgval=0.0 

return 

endif 

calculate  nominal  horizon  brightness  for  alamb 
dlamb  =alamb-280. 
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indx =int(dlamb/20.)  + 1 
indx2=indx+l 
yl=alog(hbkg(indx)) 
y2  =alog(hbkg0ndx2)) 
rati  =(amod(dlamb,20.))/20. 
yint  =  yl  +  (y2-yl)’^ratl 
hbkl  =  exp(yint) 

calculate  actual  brightness 

if(thetzr.lt.  1.5708)  then 
if(ibkflg.eq.  1)  then 

bkgval  =  hbkl  *(0.19  +  .81*sin(thetzr)) 
else 

if(ibkflg.eq.2)  then 

bkgval  =  hbkl  *(1.0  +  2*cos(thetzr)) 

else 

if(ibkflg .  eq .  3)then 

bkgval  =  1.0e-6*(hbkl*(0.19  +  .81*sin(thetzr))) 
endif 
endif 
endif 
else 

if(ibkflg.eq.l.or.ibkflg.eq.2)  then 

bkgval  =  0.1*hbkl*(1.0  +  0.7*cos(thetzr)) 

else 

bkgval  =  1.0e-6*0.1*hbkl*(l.  +  0.7*cos(thetzr)) 

endif 

endif 

return 

end 
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ATTN:  Ms.  Brathwaite 
5321  Riggs  Road 
Gaithersburg,  MD  20882 

National  Security  Agency 
ATTN:  W21  (Dr.  Longbothum) 

9800  Savage  Road 

Fort  George  G.  Meade,  MD  20755-6000 
OIC-NAVSWC 

Technical  Library  (Code  E-232) 

Silver  Springs,  MD  20903-5000 

Commander 

U.S.  Army  Research  office 
ATTN:  DRXRO-GS  (Dr.  Flood) 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27009 
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Dr.  Jerry  Davis 

North  Carolina  State  University 
Department  of  Marine,  Earth,  and 
Atmospheric  Sciences 
P.O.  Box  8208 
Raleigh,  NC  27650-8208 

Commander 

U.S.  Army  CECRL 

ATTN:  CECRL-RG  (Dr.  Boyne) 

Hanover,  NH  03755-1290 

Commanding  Officer 
U.S.  Army  ARDEC 
ATTN:  SMCAR-IMI-I,  Bldg  59 
Dover,  NJ  07806-5000 

Commander 

U.S.  Army  Satellite  Comm  Agency 

ATTN:  DRCPM-SC-3 

Fort  Monmouth,  NJ  07703-5303 

Commander 

U.S.  Army  Communications-Electronics 
Center  for  EW/RSTA 
ATTN:  AMSEL-EW-MD 
Fort  Monmouth,  NJ  07703-5303 

Commander 

U.S.  Army  Communications-Electronics 
Center  for  EW/RSTA 
ATTN:  AMSEL-EW-D 
Fort  Monmouth,  NJ  07703-5303 
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Commander 

U.S.  Army  Communications-Electronics 
Center  for  EW/RSTA 
ATTN:  AMSEL-RD-EW-SP 
Fort  Monmouth,  NJ  07703-5206 

Commander 

Department  of  the  Air  Force 
OL/A  2d  Weather  Squadron  (MAC) 

Holloman  AFB,  NM  88330-5000 

PLAVE 

Kirtland  AFB,  NM  87118-6(X)8 
Director 

U.S.  Army  TRADOC  Analysis  Center 

ATTN:  ATRC-WSS-R 

White  Sands  Missile  Range,  NM  88002-5502 

Director 

U.S.  Army  White  Sands  Missile  Range 

Technical  Library  Branch 

ATTN:  STEWS-IM-rr 

White  Sands  Missile  Range,  NM  88002 

Army  Research  Laboratory 
ATTN:  AMSRL-BE  (Mr.  Veazy) 

Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM  88(X)2-5501 

Army  Research  Laboratory 
ATTN:  AMSRL-BE-A  (Mr.  Rubio) 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM  88002-5501 
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Army  Research  Laboratory 
ATTN:  AMSRL-BE-M  (Dr.  Niles) 

Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM  88002-5501 

Army  Research  Laboratory 
ATTN:  AMSRL-BE-W  (Dr.  Seagraves) 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM  88002-5501 

USAF  Rome  Laboratory  Technical 
Library,  FL2810 
Corridor  W,  STE  262,  RL/SUL 
26  Electronics  Parkway,  Bldg  106 
GriffissAFB,  NY  13441-4514 

AFMC/DOW 

Wright-Patterson  AFB,  OH  03340-5(X)0 
Commandant 

U.S.  Army  Field  Artillery  School 
ATTN:  ATSF-TSM-TA  (Mr.  Taylor) 

Fort  Sill,  OK  73503-5600 

Commander 

U.S.  Army  Field  Artillery  School 
ATTN:  ATSF-F-FD  (Mr.  Gullion) 

Fort  Sill,  OK  73503-5600 

Commander 

Naval  Air  Development  Center 
ATTN:  A1  S'alik  (Code  5012) 

Warminister,  PA  18974 


Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-M  (Mr.  Bowers) 

Dugway,  UT  84022-5000 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-L 
Dugway,  UT  84022-5000 

Defense  Technical  Information  Center 
ATTN:  DTIC-OCP 
Cameron  Station 
Alexandria,  VA  22314-6145 

Commander 
U.S.  Army  OEC 
ATTN:  CSTE-EFS 
Park  Center  TV 
4501  Ford  Ave 
Alexandria,  VA  22302-1458 

Commanding  Officer 

U.S.  Army  Foreign  Science  &  Technology  Center 
ATTN:  CM 
220  7th  Street,  NE 
Charlottesville,  VA  22901-5396 

Naval  Surface  Weapons  Center 
Code  G63 

Dahlgren,  VA  22448-50(X) 


Commander  and  Director 
U.S.  Army  Corps  of  Engineers 
Engineer  Topographies  Laboratory 
ATTN:  ETL-GS-LB 
Fort  Belvoir,  VA  22060 

U.S.  Army  Topo  Engineering  Center 

ATTN:  CETEC-ZC 

Fort  Belvoir,  VA  22060-5546 

Commander 

USATRADOC 

ATTN:  ATCD-FA 

Fort  Monroe,  VA  23651-5170 

TAC/DOWP 

Langley  AFB,  VA  23665-5524 

Commander 
Logistics  Center 
ATTN:  ATCL-CE 
Fort  Lee,  VA  23801-6000 

Science  and  Technology 
101  Research  Drive 
Hampton,  VA  23666-1340 

Commander 

U.S.  Army  Nuclear  and  Chemical  Agency 
ATTN:  MONA-ZB,  Bldg  2073 
Springfield,  VA  22150-3198 
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James  B.  Gillespie 

U.S.  Army  Research  Laboratory 

AMSRL-IS-ES 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1197 


Patti  Gillespie 

U.S.  Army  Research  Laboratory 
AMSRL-IS-ES 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


Record  copy 

U.S.  Army  Research  Laboratory 
Attn:  AMSRL-CS-AL-TA  (Ellen  Freeman) 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


Total 


187 
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